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Abstract 

This  research  effort  concentrated  on  applying  statistical  process  control  techniques 
to  the  results  of  seven  years  worth  of  spirometry  exams  of  workers  at  Wright-Patterson 
AFB  in  helping  the  occupational  health  squadron  identify  potential  work  hazard  areas. 
Each  spirometry  exam  was  classified  as  abnormal  or  normal  based  on  a  comparison  with 
established  normals  or  a  significant  loss  of  function  from  the  previous  year’s  exam  for  the 
individual.  Each  test  was  classified  into  the  work  area  of  the  individual  and  the  number  of 
abnormalities  per  population  of  the  work  area  was  examined.  This  was  accomplished 
through  the  use  of  standardized  control  charts  for  fraction  abnormal.  The  base  population 
studied  was  “out-of-control”  for  all  seven  years;  however,  when  excluding  smokers  only 
one  year  of  the  study  for  the  whole  base  was  “out-of-control”.  Several  work  areas  were 
identified  as  being  “out-of-control”  and  recommended  for  further  study  by  the 
occupational  health  squadron. 


STATISTICAL  PROCESS  CONTROL  IN  MEDICAL  SURVEILLANCE 
AN  APPLICATION  USING  SPIROMETRY 


Chapter  1  Introduction 


1.0  Background 

Wright-Patterson  Air  Force  Base  (WPAFB)  located  in  Dayton,  Ohio  employs 
approximately  24,000  people  of  which  15,000  are  civilians  (26;2).  There  are  about  1,600 
different  facilities  located  at  WPAFB  (26:3)  and  it  is  the  responsibility  of  Medical  Group’s 
Occupational  Medicine  Element  (74th  SGPMO)  to  ensure  each  of  these  occupational 
areas  is  not  harmful  to  the  people  who  work  in  them.  A  proactive  and  preventive 
approach  to  this  goal  with  respect  to  occupational  pulmonary  diseases  is  desirable. 

In  light  of  this  goal  the  74th  SGPMO  maintains  a  health  database  on  personnel 
working  at  WPAFB.  One  to  seven  years  worth  of  pulmonary  function  data  exists  on  each 
of  the  individuals  monitored;  the  database  also  includes  their  work  area  (building  and 
office),  a  family  history  of  previous  health  problems,  and  personal  habits  of  the  individuals 
(i.e.  cigarette  smoking  and  alcohol  consumption).  However,  the  74th  SGPMO  only 
maintains  this  database  and  does  not  conduct  analysis.  The  purpose  of  this  research  is  to 
conduct  the  needed  analysis  to  help  the  74th  SGPMO  monitor  the  different  work  zones. 

The  data  extracted  from  the  database  for  this  analysis  includes  the  following  as 
identifiers:  Social  Security  Number  and  Work  Zone  history.  Independent  variables 
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include:  Date  of  Birth,  Gender,  Race,  Height,  family  history  information  on  Blood 
Disease,  Asthma,  Lung  Disease,  Liver  Disease,  Hepatitis,  and  Jaundice,  personal 
information  on  packs  of  cigarettes  smoked  per  day,  number  of  years  smoked  if  quit,  and 
total  number  of  years  smoked  if  currently  smoking.  Also  ounces  of  liquor,  bottles  of  beer, 
and  number  of  wine  glasses  consumed  per  week  is  included.  Possible  dependent  variables 
are  the  results  of  the  pulmonary  function  test;  these  include  Forced  Vital  Capacity 
(FVC),  Forced  Expiratory  Volume  in  One  Second  (FEVi),  the  Ratio  of  FEVi  /  FVC,  the 
percent  of  predicted  FVC  (FVCpred),  the  percent  of  predicted  FEVi  (FEVlpred),  and  the 
percent  of  predicted  Ratio  FEVi  /  FVC  (Ratpred).  The  number  of  individuals  who  have 
lung  data  monitored  is  1,945  (as  of  3  October  1995).  The  following  table  shows  the 
breakdown  on  number  of  measurements  per  person  in  the  database. 


Table  1-1 

Pulmonary  Function  Tests  Measurements 


Number  of  Pulmonary  Function  Tests 

Number  of  People 

1 

815 

2 

469 

3 

248 

4 

172 

5 

104 

6 

110 

7 

27 

The  pulmonary  function  test  is  accomplished  on  a  spirometer.  A  spirometer  is  “an 
instrument  for  measuring  air  entering  and  leaving  the  lungs”  (25;  1 122).  The  basic 
principle  is  that  the  subject  inhales  as  much  air  as  possible  and  then  forcefully  expires  the 
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air  as  hard  as  possible  as  long  as  possible  into  the  spirometer.  The  spirometer  measures 
the  respective  variables  during  the  process.  The  standard  used  in  this  case  is  three  trials 
and  the  “best  effort”  is  recorded  as  the  final  result.  The  type  of  spirometer  used  in  this 
analysis  is  a  SPIROTECH  S400  VER  6C  Andersen/Spirotech  manufactured  by  Spirotech, 
a  division  of  Graseby  Andersen..  Figure  1-1  shows  a  type  of  spirometer  similar  to  the  one 
used  in  this  study.  According  to  Mr.  William  M.  Yancey,  the  Pulmonary  Function  and 
ECG  Technician  at  the  74th  SGPMO,  it  is  approved  and  does  meet  the  standards  set  for 
spirometers  set  by  the  American  Thoracic  Society  (ATS).  Please  see  Figure  1-2  for  the 
sample  output  of  such  a  test.  As  stated  above,  the  measurements  used  in  this  study  are 
FVC,  FEVi,  FEVi  /  FVC  and  the  corresponding  predicted  values.  FVC  is  the  maximum 
amount  of  air  which  can  be  exhaled  with  a  maximum  forced  effort  from  a  position  of 
maximum  inspiration  (12:17  &  15;A4).  It  is  the  most  air  a  person  can  breath  out  after 
inhaling  as  much  as  possible.  FEVi  is  the  volume  of  air  exhaled  in  the  first  second  of  the 
FVC  performance  (12:18  &  15:A7).  Both  FVC  and  FEVi  are  measured  in  liters  BTPS, 
Body  conditions:  normal  body  temperature;  37  degrees  Celsius,  ambient  pressure, 
saturated  with  water  vapor  (2:2-3).  The  ratio  is  simply  the  volume  of  air  exhaled  in  the 
first  second  divided  by  the  FVC  value  and  converted  to  a  percentage  by  multiplying  by 
100.  The  predicted  values  of  FVC,  FEVi,  and  FEVi  /  FVC  are  calculated  by  the 
spirometer  based  on  Knudson  normals  (which  will  be  defined  in  Chapter  2). 

Normal  values  and  variations  among  the  spirometry  results  for  FVC  are  based  on 
Gender,  Age,  Height,  Weight,  Race,  Technical  (due  to  the  machine)  and  Unexplained 
factors  (1 :2  &  7:325).  Please  see  Table  1-2.  Studies  show  black  adults  values  of  FVC 
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Figure  1-1 

A  Spirotech  Spirometer 


and  FEVi  range  from  5  to  20%  lower  than  White  adults  (5:251).  Blacks  of  the  same 
height  as  Whites,  in  general,  have  longer  legs  and  a  shorter  thorax;  therefore  these  values 
tend  to  be  lower  when  controlling  for  height  (10:828).  The  ratio  of  FEVi  /  FVC,  is  also 
higher  in  non-whites  because  the  reduction  in  FEVi  is  less  than  the  reduction  in  FVC  for 
these  subjects  (5:250).  Technical  factors  include  the  instrument  being  used,  the 
performance  of  the  test  administrator,  the  performance  of  the  subject,  interaction  between 
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the  administrator  and  subject,  the  procedure  being  used,  and  the  temperature  /  altitude  the 
test  was  administered  at  (5:244).  The  unexplained  factors  include  past  and  present  health 
experiences  (personal  health  history,  family  health  history,  and  current  health),  past  and 
present  exposures  (i.e.  cigarette  smoking  and  occupational  exposures)  and  socioeconomic 
factors  (7:325). 


Table  1-2 

The  Proportion  of  Variance  in  FVC,  or  FEVi 
(1:2  &  7:325) 


Factor 

ProDortion  of  Exolained  Variance 

Gender 

0.30 

Age 

0.08 

Height 

0.20 

Weight 

0.02 

Ethnic  Differences 

0.10 

Technical 

0.03 

Unexplained 

0.27 

1 . 1  Statement  of  the  Problem 

Using  the  information  extracted  from  the  health  database,  accomplish  an  analysis  to 
determine  any  signs  of  possible  adverse  effects  to  individual’s  pulmonary  functions  which 
may  be  a  result  of  their  work  environment.  The  scope  of  this  effort  is  to  determine  if  there 
is  a  possible  adverse  effect  between  a  work  zone  and  the  employee’s  health  based  on  the 
spirometry  results.  This  analysis  is  to  be  used  as  a  screening  tool  for  the  74th  SGPMO.  If 
a  work  zone  does  not  meet  “normal”  criteria,  it  is  the  responsibility  of  the  74th  SGPMO  to 
determine  if  corrective  action  is  warranted  for  that  zone. 
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SPIROTECH  S400  VER  6C 

Andersen/Spirotech 

Atlanta,  Georgia 

ID:  123456789 

Time:  9:38 

Date:  6/  8/95 

Age:  27  yrs 

Height:  72  ins 

Sex:  Male 

Race: 

WHITE 

22.  rc 

760.0  mmHg  BTPS 

:  1.090 

Normals:  Knudson  Chemiack 

Best  Effort:  Test  #3 

Criteria  met 

Parameter 

Actual 

Predicted 

%Predicted 

FVC 

5.66 

5.65 

100% 

FEV.5 

3.33 

3.56 

93% 

FEVl 

4.40 

4.58 

96% 

FEV3 

5.47 

5.44 

101% 

FEF25-75 

3.82 

5.53 

69% 

FEF75-85 

1.27 

PEFR 

10.56 

10.25 

103% 

FEF25 

7.35 

9.53 

77% 

FEF50 

4.33 

6.81 

63% 

FEF75 

1.82 

3.58 

51% 

FEV.5/FVC 

.59 

FEVl/FVC 

.78 

.84 

93% 

FEV3/FVC 

.97 

FET 

7.84 

Figure  1-2 

Sample  Output  of  Spirometry  Test 
The  Units  used  are  liters  BTPS 
The  Predicted  Column  is  the  Knudson  Normal 


Chapter  2  Literature  Review 


2, 1  Spirometry  and  Occupational  Health 

This  section  will  deal  with  the  justification  for  using  the  spirometry  results  for 
analysis  of  possible  occupational  cause  of  pulmonary  disease.  The  Occupational  Health 
and  Safety  Administration  (OSHA)  requires  spirometry  for  people  working  in  certain 
exposure  environments  (12:2,  13:229).  These  include  employee  exposure  to  asbestos, 
coke  oven  emissions,  and  cotton  dust  (14:229).  The  National  Institute  for  Occupational 
Safety  and  Health  (NIOSH)  recommends  pulmonary  function  testing  for  exposure  to 
numerous  other  substances;  beryllium,  cadmium,  formaldehyde,  nitrogen-oxides,  etc... 
(14:229).  In  Albert  Miller’s  “Application  of  Pulmonary  Function  Tests”  article,  he  states 
pulmonary  function  tests  may  be  the  only  evidence  of  certain  types  of  diseases  and  in  the 
case  of  obstructive  related  diseases,  they  are  actually  better  than  a  chest  x-ray  (16:4). 

Even  though  the  results  of  the  spirometry  tests  do  not  establish  the  diagnosis,  they  are 
helpful  in  suggesting  a  possible  abnormality  (16:4).  In  McKay’s  article  on  “Pulmonary 
Function  Testing  in  Industry”,  he  states  spirometry  is  now  regarded  as  a  key  part  of  any 
respiratory  surveillance  program  (14:229).  Not  only  is  spirometry  a  good  instrument  with 
respect  to  screening  for  lung  disease  (4:349),  but  it  can  identify  pulmonary  function 
abnormalities  which  may  be  overlooked  otherwise  (6:25).  Now,  why  of  all  the  spirometry 
output  shown  in  Figure  1-1,  only  three  variables  are  to  be  used  in  this  study? 

In  his  article  “Respiratory  Disorders”,  Dr.  Wegman  states  the  most  useful 
information  from  the  spirometry  test  is  the  results  of  the  FVC,  FEVi,  and  the  FEVi  /  FVC 
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ratio  when  evaluating  possible  occupational-related  respiratory  illnesses  (24:321).  In  the 
same  article  he  also  provides  a  table  of  major  occupational  pulmonary  diseases  and  how 
these  three  variables  react  in  each  of  the  situations.  Please  see  Table  2-1  for  a  summary  of 
the  different  diseases. 


Table  2-1  (24:326) 

Impact  of  Disease  Type  on  Variables 


Type 

Occupational  Example 

FEVl 

FVC 

Ratio  FEVl/FVC 

Restrictive 

Silicosis, 

Asbestosis 

Decrease 

Decrease 

No  Change 

Obstructive 

Byssinosis, 

isocyanate  asthma 

Decrease 

No 

Change 

Decrease 

Granulomata 

Beryllium  disease 

Decrease 

Decrease 

N/A 

Pulmonary 

edema 

Cadmium  Poisoning 

Decrease 

Decrease 

N/A 

Restrictive  lung  disease  means  the  lungs  are  too  stiff  and/or  too  small.  An 
individual  with  restrictive  lung  disease  will  slowly  experience  difficulty  in  breathing 
(9: 175).  An  obstructive  lung  disease  (a.k.a.  Chronic  Obstructive  Pulmonary  Disease) 
means  the  airways  are  becoming  blocked  and/or  the  lungs  are  losing  their  structural 
integrity.  With  respect  to  occupational  asthma,  the  airways  become  inflamed  due  to  the 
irritant  and  this  results  in  blocking  (9: 145).  Granulomata  indicates  the  presence  of  lung 
nodules.  A  pulmonary  edema  is  basically  fluid  in  the  lungs. 
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The  ATS  reinforces  Dr.  Wegman’s  view  by  stating  the  FVC,  FEVi,  and  FEVj  / 
FVC  ratios  are  the  three  basic  factors  in  the  interpretation  of  a  spirometry  exam  and  goes 
on  to  state  in  diagnosing  for  an  obstructive  disease,  the  FEVi  /  FVC  ratio  is  the  most 
important  measurement  (1 : 1212).  Several  studies  on  occupational  lung  disease  utilized 
these  three  variables;  therefore  setting  precedence  for  the  approach  taken  in  this  research. 
Parker’s  study  on  the  pulmonary  function  of  autobody  repair  workers  was  based  on 
measurements  of  FEVi,  FVC  and  the  FEVi  /  FVC  ratio;  the  FEVi  /  FVC  ratio  results 
were  used  to  show  indications  of  possible  obstructive  disease  (22:768-771).  Monson  cites 
a  study  of  Potash  mine  workers  in  Canada  in  which  the  pulmonary  function  was  measured 
by  spirometry.  He  showed  how  the  FEVi  value  related  to  SO2  exposure  (19:215-216). 

2  2  Defining  “Abnormality”  With  Respect  to  a  Reference  Population 

As  stated  in  Chapter  1,  general  population  normal  values  for  FVC,  FEVi,  and  the 
FEVi  /  FVC  ratio  exist  with  respect  to  the  factors  sex,  age,  height,  weight,  and  race.  The 
analysis  of  a  specific  population  is  possible  because  of  these  general  population  normal 
values  (19:214).  The  spirometer  used  in  this  study  outputs  a  percent  of  predicted  value 
for  each  person  based  on  the  Knudson  method.  Table  2-2  shows  the  coefficients  for  the 
regression  equations  the  Spirotech  spirometer  uses  in  predicting  the  normal  values  based 
on  the  equation: 

y  =  C  +  [age(years)  *  age  coefficient]  +  [standing  ht(cm)  *  ht  coefficient]  (13:589) 
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where  y  is  either  the  predicted  FVC  value,  the  predicted  FEVi  value,  or  the  predicted 
FEVi  /  FVC  ratio  value.  For  blacks,  the  predicted  values  for  FVC  and  FEVi  need  to  be 
multiplied  by  0.85  to  account  for  the  racial  difference  (10:848). 

In  the  article  “Predictable  Confusion”,  Dr.  Glindmeyer  states  one  set  of  regression 
equations  could  show  a  person  to  be  normal  with  respect  to  their  lung  function,  while 
another  set  would  show  the  person  to  be  abnormal  (10:845).  He  recommends  using  the 
values  of  Knudson  because  they  are  based  on  a  large  sample  of  healthy,  never-smokers, 
who  are  free  from  any  respiratory  symptoms  (10:848).  This  is  consistent  with  the 
recommendations  of  OSHA  which  adopted  the  predicted  normals  based  on  Knudson’ s 
testing  for  its  1978  cotton  dust  standard  (12:2). 


Table  2-2 

Coefficients  for  Knudson  Regression  Equations  (13:590) 


Test 

Gender 

Age 

C 

Coefficient 

Age 

Coefficient 

Ht 

Coefficient 

Std  Dev  of 

Prediction 

FVC 

Male 

<  25  Years 

-5.508 

0.078 

0.050 

0.544 

FVC 

Male 

>  25  Years 

-5.459 

-0.029 

0.065 

0.601 

FVC 

Female 

<  20  Years 

-3.469 

0.092 

0.033 

0.500 

FVC 

Female 

>  20  Years 

-1.774 

-0.022 

0.037 

0.519 

FEVi 

Male 

<  25  Years 

-4.808 

0.045 

0.046 

0.523 

FEVi 

Male 

>  25  Years 

-4.203 

-0.027 

0.052 

0.541 

FEVi 

Female 

<  25  Years 

-2.703 

0.085 

0.027 

0.422 

FEVi 

Female 

>  25  Years 

-0.794 

-0.021 

0.027 

0.434 

FEVi/FVC 

Male 

All  Ages 

103.64 

-0.140 

-0.087 

6.721 

FEVi/FVC 

Female 

All  Ages 

107.38 

-0.109 

-0.111 

7.664 
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In  Bascom  and  Ford’s  article,  “Don’t  Just  ‘Do  Spirometry’-  Closing  the  Loop  in 
Workplace  Spirometry  Programs”,  the  answer  to  the  question,  “Is  the  output  normal?”,  is 
answered  by  the  comparison  of  the  individual  test  result  with  the  predicted  value  based  on 
the  regression  equation  used  (4:355).  The  general  cutoff  for  abnormals  recommended  for 
FEVi  and  FVC  is  less  than  80%  of  the  predicted  value  (4:355).  The  standard  deviations 
about  the  regression  line  for  these  predictions  are  shown  in  Table  2-2.  In  Miller’s 
“Prediction  Equations  and  ‘Normal’”  article,  he  cites  the  American  Lung  Association’s 
handbook  on  Chronic  Obstructive  Pulmonary  Disease  recommendation  that  80%  of 
predicted  is  abnormal  with  respect  to  both  the  FEVi  and  FVC  variables  (17:203).  With 
the  prediction  equations,  several  organizations  set  guidelines  for  determining  abnormality 
of  the  pulmonary  function  with  respect  to  the  predicted  values.  The  NIOSH,  ATS,  and 
American  Medical  Association  (AMA)  recommend  the  lower  limits  of  Normal  as  shown  in 
Table  2-3. 


Table  2-3  (15:C-7) 

Recommended  Lower  Limits  of  Normal 


FVC  (pet  of 

predicted  value) 

FEVI  (pet  of 

predicted  value) 

FEVI /FVC 

Observed  Value 

NIOSH 

80 

80 

70 

AMA 

80 

80 

70 

ATS 

80 

80 

75 

11 


The  American  Thoracic  Society  does  recommend  another  approach  to  defining  an 
abnormal  reading  with  respect  to  the  population  predicted  values.  They  recommend 
normal  ranges  based  on  fifth  percentiles  of  the  reference  population,  which  is  considered 
to  be  below  the  lower  limit  of  normal  (1 .1206)  Miller  agrees  with  this  practice  because 
having  a  fixed  percent  predicted  value  will  cause  false  negatives;  classification  of  “normal” 
readings  as  “abnormal”  (17205).  Dockery  adds  that  the  use  of  80%  of  the  predicted  value 
as  a  lower  limit  will  more  likely  misclassify  shorter  and  older  people  as  “abnormal”  and 
also  misclassify  taller  and  younger  people  as  “normal”  (7;326).  This  is  because  the 
variation  by  height  and  age  is  not  consistent.  McKay  acknowledges  the  80%  lower  limit 
for  normal  value,  but  believes  a  95th  percentile  method  is  more  valid  (14:232)  because  this 
method  will  have  fewer  false  positives.  Table  2-4  summarizes  the  abnormality  cutoff  (the 
lower  fifth  percentile)  with  respect  to  the  general  population  normal  values  using  the 
Knudson  regression  equations  and  the  percent  of  predicted  values. 


Table  2-4  (15:C-11,12) 

Lower  Fifth  Percentiles  Based  on  Knudson 
Values  are  Percent  of  Predicted 


Gender 

Age 

FVC 

FEV, 

FEV.  /  FVC  Ratio 

Male 

25-39  years 

81.1% 

79.1% 

87.0% 

40-84  years 

73.4% 

77.2% 

87.0% 

Female 

20-39  years 

76.9% 

70.3% 

85.4% 

40-70  years 

75.2% 

77.9% 

85.4% 
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2.3  Defining  Abnormal  Declines  of  Lung  Function 


The  database  does  have  multiple  readings  of  individuals,  therefore  a  longitudinal 
study  can  also  be  accomplished  with  respect  to  this  data.  According  to  Dr.  McKay,  even 
though  we  can  compare  a  person’s  test  results  to  the  reference  values,  comparing  them  to 
previous  tests  is  more  desirable  (14:234).  This  is  because  the  coefficient  of  variation  of  a 
given  test  within  the  individual  is  smaller  than  the  population’s  coefficient  of  variation 
(14:234).  The  coefficient  of  variation  is  the  standard  deviation  as  a  proportion  of  the 
mean  of  the  population.  Where  the  population’s  variation  is  based  on  all  the  factors  listed 
in  Table  1-2,  the  individual’s  variation  will  be  limited  to  the  technical  and  unexplained 
factors  of  taking  the  spirometry  exam  itself  Bascom  and  Ford  state  using  the  person’s 
own  longitudinal  data  is  much  better  than  comparing  a  single  value  versus  a  population 
predicted  value  (4:357).  In  the  event  of  having  multiple  test  results  letting  the  person  be 
their  own  control  is  much  more  desirable  (12:29)  Table  2-4  shows  the  percent  changes 
before  a  meaningful  diagnosis  of  abnormality  for  a  longitudinal  study  can  be  claimed. 


Table  2-4  (15:H-6,7) 
Abnormality  in  a  Longitudinal  Study 


Source 

FVC 

FEV, 

FEV, /FVC Ratio 

ATS 

Annual  Decline  of 

Annual  Decline  of 

15%  or  more 

15%  or  more 

NIOSH 

Annual  Decline  of 

Annual  Decline  of 

Annual  Decline  of 

10%  or  more 

10%  or  more 

5%  or  more 
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2.4  Summary  of  Classification 


From  the  information  stated  above,  the  cuts  for  an  “abnormality”  classification 
with  respect  to  the  result  of  the  spirometry  test  will  be  a  combination  of  the  reference 
population  normals  using  5%  of  the  lower  limit  (based  on  the  values  in  Table  2-3)  for 
FVC  and  FEVi  and  a  longitudinal  study  where  multiple  tests  on  individuals  is  available 
(based  on  the  values  in  Table  2-4).  For  ratio  the  observed  value  of  0.70  or  lower  will  be 
used  and  a  longitudinal  decline  of  10%  or  more.  The  5%  cutoff  yielded  too  many  false 
abnormalities.  If  the  person  falls  into  any  of  these  categories  their  classification  will  be 
abnormal.  Dr.  Glindmeyer  states  a  subject  should  be  classified  as  abnormal  if  he/she  is 
below  the  normal  range  or  exhibits  an  abnormal  longitudinal  decline,  even  if  it  is  within  the 
normal  range  (12:849).  Another  reason  for  combining  the  two  methods  in  this  study  is  the 
information  is  available  and  if  the  person  enters  the  work  environment  with  a  well  above 
predicted  value  and  declines  at  20%  per  year,  a  study  with  respect  to  only  the  reference 
population  values  would  not  classify  this  person  as  abnormal  until  it  is  too  late  for  the 
screening  purpose  of  the  study  to  be  useful  (4:357). 

2.5  Cigarette  Smoking 

No  epidemiological  study  on  pulmonary  disease  can  be  accomplished  without 
accounting  for  cigarette  smoking  (1:1206,  19:214).  It  is  definitely  a  potential  confounding 
factor  in  the  analysis  of  possible  pulmonary  disease  for  a  work  group  area  (1 : 1206). 
Cigarette  smoking  does  lead  to  airflow  obstruction,  which  will  show  in  a  lower  FEVi 
value  (6:25).  Even  if  a  smoker  has  one  “normal”  pulmonary  function  test,  he/she  should 
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not  believe  they  are  free  of  problems;  their  loss  of  lung  function  will  show  in  the  future 
(6;25).  Almost  all  of  the  literature  states  smoking  does  need  to  be  accounted  for  in  this 
study.  The  ATS  states  it  does  need  to  be  an  independent  variable  in  the  analysis,  but  the 
most  appropriate  method  of  measurement  is  not  known ;  i.e.  binary  yes/no  variable, 
amount  currently  smoked,  pack-years,  etc.(l :  1206).  The  three  variables  in  the  database 
are:  current  smoking  status,  past  smoking  status,  and  number  of  packs  smoked  per  day. 

2.6  Statistical  Process  Control 

Each  observation  of  an  individual  will  be  classified  as  “normal”  or  “abnormal”  with 
respect  to  the  FVC,  FEVi,  and  the  ratio  FEVi  /  FVC.  A  table  showing  the  logic  for  the 
classifications  is  presented  in  chapter  3.  These  classifications  are  the  attributes  of  the 
analysis.  The  statistical  quality  technique  which  corresponds  to  the  fraction  of  abnormal 
occurrences  per  work  zone  is  a  control  chart  for  fraction  nonconforming;  a  p  chart 
(20: 147).  According  to  Montgomery,  the  item  may  have  several  quality  characteristics  (in 
this  case  each  of  the  individual  FVC,  FEVi,  and  ratio  FEVi  /  FVC  readings  as  well  as  the 
longitudinal  evaluations)  which  may  be  classified  as  abnormal.  If  any  of  these  are 
abnormal,  the  subject  is  classified  as  abnormal  (20:148).  The  fraction  defined  as  abnormal 
would  then  be  the  ratio  of  the  number  of  people  in  each  work  area  classified  as  abnormal 
divided  by  the  total  number  of  people  in  that  work  area  (20: 148).  Based  on  this  fraction,  a 
standard  value  of p  decided  by  management  (the  74th  SGPMO)  can  be  used  to  develop 
the  upper  control  limit  (UCL),  the  center  line,  and  the  lower  control  limit  (LCL)  in  a 
control  chart  (20: 149).  The  UCL  and  the  LCL  define  the  range  in  which  the  fraction 


15 


classified  as  abnormal  would  be  accepted,  or  in  other  words,  the  work  area  would  be 
considered  in  control  with  respect  to  lung  abnormalities.  The  equations  for  the  UCL, 
center  line,  and  LCL  follow: 


C/CL-p  +  3  1^^^ 

(20:149) 

V  w 

Center  line  =  p 

(20:149) 

LCL  =  p 

V  « 

(20:149) 

In  the  above  equations  ‘p’  is  the  fraction  nonconforming  decided  by  management 
and  ‘n’  is  the  number  of  individuals  tested  in  the  work  area.  When  the  fraction  of 
abnormalities  is  plotted  for  each  work  area,  the  points  which  fall  above  the  UCL  identify 
the  work  areas  for  further  investigation  to  determine  if  there  is  an  assignable  cause 
(20:153).  In  the  case  where  sample  sizes  vary  (which  they  do  for  the  number  of  people 
per  work  area)  the  simple  approach  is  to  base  the  control  limits  on  each  respective  sample 
size  (20: 163).  This  is  the  method  employed  in  this  study.  For  presentation  purposes  the 
“standardized”  control  chart  is  used.  This  plots  the  same  information  as  the  p  chart; 
however,  in  standard  deviation  units  (20: 167).  The  UCL  is  3  standard  deviation  units  and 
the  variable  which  is  plotted  is  Z  and  is  found  by  the  following  equation: 


Z  = 


ip-p) 

\P(}-P) 


(20:149) 


V  n 

where  ‘  Z  ’  is  the  Sample  Standard  Deviation,  ‘  ^  ’  is  the  sample  percent  abnormal,  and  ‘p’ 
is  the  Center  Line  from  above. 
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Chapter  3  Methodology 


3.1  Overview 

The  steps  in  transferring  the  database  from  the  Phoenix  system  to  the  SAS  system, 
the  logic  used  in  creating  a  workable  database,  and  the  programs  used  in  the  analysis  will 
be  outlined  in  this  chapter.  This  thesis  effort  was  accomplished  in  conjunction  with  the 
same  type  of  research  with  respect  to  liver  tests,  therefore  the  logic  for  creating  the 
database  including  the  liver  data  is  also  presented.  Figure  3-1  depicts  the  process  flow  in 
the  creation  of  a  SAS  database,  the  classification  of  abnormalities,  and  the  implementation 
of  control  charts  in  which  the  lung  and  work  zone  data  can  be  analyzed. 


3.2  Creating  the  SAS  Database 

The  first  step  involves  downloading  the  appropriate  data  elements  from  the 

Phoenix  Health  database  system  to  flat  ASCII  files.  This  was  accomplished  using  the 

Sentry  Health  Surveillance  system  (a  menu  driven  data  inquiry  system  for  the  Phoenix 

database).  Due  to  the  structure  of  this  database  the  following  7  output  files  were  created: 

Thesisl:(Information  on  Date  of  Birth,  Sex,  Race,  Height,  Date  of  Exam) 

ssan _ dob  sex  race  height  examdate 

007645826  1961  M  C  72  1995/05/26 

Thesis2:(Lung  Data  Information) 


ssan 

date 

fevl 

fvc 

ratio 

fvcpct  fevpct  ratiop 

002547331 

1993/11/17 

3.01 

3.61 

0.83 

98  100  98 

Thesis3:(Liver  Data  Information) 


sgot 

date 

ap 

bilirubin 

albumin 

005527773  25 

27 

1989/11/06 

96 

106 

0.5 

4.3 
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Figure  3-1 
Methodology 


Thesis4:(White  Blood  Count  Data  Information) 

ssan _ wbc  date  wbc  hematocrit 

002547331  7.2  1993/05/13  38.9 

Thesis5:(Work  Zone  Information) 

ssan _ zone _ start  date  end  date 

002485546  B640D3  1993/02/08  1994/04/12 

Thesis6:(Family  History  Information  on  Blood  Disease,  Asthma,  Lung  Disease, 
Liver  Disease,  Hepatitis,  and  Jaundice;  Personal  Information  on  Packs  of  day 
smoked.  Number  of  years,  and  how  long  since  they  quit) 


ssan  bid  ast 

lun 

liv 

hep 

iau 

smo 

smo 

smo 

005208011  B  B 

B 

B 

B 

B 

2 

20 

30 

A  “B”  entry  means  no  history  and  an  “A”  entry  means  there  is  a  history. 

Thesis7:(Personal  information  on  ounces  of  liquor,  bottles  of  beer,  and  number  of 
wine  glasses  consumed  per  week) 

ssan _ dl  d2  d3 

033420318  223 
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As  each  exam  is  repeated,  a  new  entry  based  on  the  SSAN  is  created;  therefore  we 
do  have  multiple  entries  on  the  same  person.  After  the  downloading  of  each  respective  file 
to  floppy  disks,  the  files  are  transferred  to  the  UNIX  mainframe  system  at  AFIT  via  the 
WS-FTP  protocol. 

The  third  step  in  this  process  is  creating  a  SAS  database  of  the  above  information 
in  a  workable  format.  There  are  several  steps  involved  in  this  process.  The  first  one  is 
making  a  SAS  database  out  of  the  above  files.  This  is  accomplished  in  the 
CONVERT.  SAS  program  (please  see  appendix  I).  After  each  of  the  databases  are 
created,  they  are  then  converted  back  to  flat  ASCII  files  so  that  the  information  (FVC, 
FEVi,  SGPT,  etc. )  is  on  the  same  line  as  the  SSAN.  In  this  step,  the  multiple 
observations  per  individual  is  done  away  with  and  there  now  exists  flat  files  where  all  the 
information  per  SSAN  is  on  one  line.  This  is  accomplished  in  the  *RAW.SAS  programs 
(please  see  appendix  II).  After  each  *RAW.SAS  program  runs,  the  output  file  is  saved  as 
a  flat  ASCII  file.  The  next  step  involves  reading  back  in  the  raw  files  into  SAS  databases 
and  making  one  dataset  for  the  entire  information.  This  is  accomplished  in  the 
MERGEALL.SAS  program  (please  see  appendix  III).  This  creates  a  database  with  all  the 
lung,  liver,  and  zone  information  which  is  called  “HEALTH. WPAFB2”. 

3,3  Classification  of  Common  Work  Zone  Exposure  Areas 

The  work  zones  are  based  on  the  area  of  base  in  which  a  person  works;  this  can  be 
either  A,  B,  C,  or  K,  the  building  number  in  which  they  work,  a  letter  for  identifying 
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common  exposures,  and  a  number  for  furthur  breakdown  of  the  exact  area  within  the 
building  they  work.  For  purposes  of  this  analysis,  common  exposure  areas  for  work  zones 
will  be  based  on  the  area,  building  number,  and  first  letter  of  exposure. 


3.4  Classification  of  Abnormal  Lung  Function  Test 

Step  four  involves  classifying  each  of  the  individual  FVC,  FEVi,  and  FEVi  /  FVC 
readings  as  either  normal  or  abnormal.  The  LUNGYRSl.SAS  and  LUNGYRS2.SAS 
programs  (please  see  appendix  IV)  accomplish  this  task  while  also  ensuring  the  correct 
classification  of  zone  for  each  observation  is  made.  The  logic  for  classifying  a  test  as 
abnormal  based  on  the  literature  review  in  chapter  2  is  summarized  in  Table  3-1. 

Table  3-1 


Logic  Table  for  Classification  of  Abnormalities 


Test 

Gender 

Age 

Value 

Pet  Predicted  of  FVC 

Male 

<  40  years 

81.1  or  less 

Pet  Predicted  of  FVC 

Male 

>  39  years 

73.4  or  less 

Pet  Predicted  of  FVC 

Female 

<  40  years 

' 

76.9  or  less 

Pet  Predicted  of  FVC 

Female 

>  39  years 

75.2  or  less 

Pet  Predicted  of  FEVi 

Male 

<  40  years 

79.1  or  less 

Pet  Predicted  of  FEVi 

Male 

>39  years 

77.2  or  less 

Pet  Predicted  of  FEVi 

Female 

<  40  years 

70.3  or  less 

Pet  Predicted  of  FEVi 

Female 

>  39  years 

77.9  or  less 

Observed  Ratio 

Both 

All  Ages 

70  or  less 

Annual  Decline  of  FVC 

Both 

All  Ages 

1 5%  or  more 

Annual  Decline  of  FEVi 

Both 

All  Ages 

15%  or  more 

Annual  Decline  of  Ratio 

Both 

All  Ages 

10%  or  more 
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The  following  assumptions  and  deletions  were  made: 

-  Only  one  pulmonary  function  measurement  per  year  per  SSAN  is  allowed 

-  If  the  date  of  the  first  pulmonary  function  test  is  missing  and  a  second  test  exists, 
it  is  made  exactly  one  year  before  the  date  of  the  second  test. 

-  If  date  of  birth  was  a  missing  value,  the  subject  was  classified  as  over  40  years 

-  If  gender  was  missing,  the  subject  was  classified  as  a  Male 

-  If  the  subject  had  multiple  lung  exams  in  one  year,  only  the  last  exam  is  counted 

3,5  Creating  Control  Charts 

The  output  data  fi-om  the  LUNGYRS2.SAS  program  will  be  saved  to  a  raw  file 
and  imported  into  a  Microsoft  Excel  Spreadsheet.  The  Spreadsheet  database  has  the 
format  displayed  in  Table  3-2. 

After  the  data  is  input  into  spreadsheet  format, /^-charts  and  standardized  control 
charts  will  be  developed  using  the  formulas  in  chapter  2.  Two  main  charts  will  be  created. 
The  first  one  will  be  a  p-chdxi  with  a  center  line  of  0. 10  measuring  the  percent  of 
abnormalities  for  all  people  in  all  zones.  The  second  chart  will  also  be  a  /?-chart  with  a 
center  line  of  0. 10.  This  chart  will  exclude  all  people  who  smoke  and  will  only  measure 
abnormality  rates  of  nonsmokers  in  each  respective  zone.  A  column  included  in  the 
spreadsheet  will  identify  the  case  numbers  which  fall  above  the  Upper  Control  Limit 
indicating  an  “out-of-control”  process.  This  information  will  then  be  summarized  and 
given  to  the  74th  SGPMO  for  action. 
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Table  3-2 

Excel  Spreadsheet  Database 


Data  Name 

Description 

CASE 

Zone  Identifier  by  Case  Numbers 

ABNML2 

Total  Number  of  Nonsmoker  Abnormalities  in  Zone 

ABNORMAL 

Total  Number  of  Abnormalities  in  Zone 

COUOT 

Total  Number  of  People  in  Zone 

COUNT2 

Total  Number  of  Nonsmokers  in  Zone 

SAB 

Total  Number  of  Smoker  Abnormalities  in  Zone 

SMOKE 

Total  Number  of  Smokers  in  Zone 

YEAR 

Year  of  Test 

ZONE 

Work  Area  Identifier 

CONTROL 

1  if  out  of  control,  0  if  in  control 
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Chapter  4  Results 


4.1  Overview  of  Chapter 

The  results  of  this  analysis  effort  will  be  summarized  in  this  chapter.  The  first  part 
will  include  a  summary  of  the  number  and  type  of  abnormalities  for  each  of  the  years  1989 
through  1995.  The  second  part  will  include  standardized  control  charts  and  /7-charts  to 
graphically  indicate  which  work  areas  and  years  are  considered  above  an  acceptable  limit 
with  respect  to  the  number  of  abnormal  lung  results.  This  section  will  also  include  a 
summary  of  the  work  areas  which  were  consistently  out  of  control. 

4.2  Summary  of  Abnormal  Classifications 

There  exists  up  to  six  possible  ways  for  an  individual’s  lung  function  test  to  be 
classified  as  abnormal  using  the  logic  in  Table  3-1.  The  FVC  or  FEVi  may  be  less  than  a 
certain  percentage  of  a  predicted  value,  the  FEVi  /  FVC  ratio  may  be  less  than  a  certain 
percentage,  or  any  of  the  three  measurements  may  have  a  larger  than  expected  normal 
decline  fi-om  the  previous  year’s  measurement.  Table  4-1  shows  the  number  of  people 
with  at  least  one  abnormal  reading  from  the  six  possibilities,  the  total  number  of  individual 
tests’  administered  in  each  year,  the  number  of  nonsmokers,  the  number  of  nonsmokers 
with  at  least  one  abnormal  reading  and  the  percentage  of  abnormal  readings  by  smokers 
and  nonsmokers  for  each  year  of  the  study.  An  interesting  fact  is  that  the  percentage  of 
nonsmokers  with  an  abnormal  lung  function  test  is  lower  than  the  population  group  as  a 
whole  for  each  of  the  years  under  study. 
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Table  4-1 

Summary  of  Abnormal  Tests  by  Year 


Year 

Total 

Number 

Abnormal 

Pet 

Abnormal 

Total 

Nonsmokers 

Number 

Abnormal 

Nonsmokers 

Pet 

Abnormal 

Nonsmokers 

1989 

232 

59 

25% 

94 

18 

19% 

1990 

923 

215 

23% 

443 

96 

22% 

1991 

854 

137 

16% 

420 

50 

12% 

1992 

733 

105 

14% 

311 

31 

10% 

1993 

836 

127 

15% 

379 

45 

12% 

1994 

676 

104 

15% 

316 

31 

10% 

1995 

290 

54 

19% 

127 

16 

13% 

Table  4-2  shows  the  frequency  of  the  abnormal  “hits”  for  an  individual  lung  test 
for  each  year  of  the  study.  The  range  of  values  is  from  one  to  six  as  described  above. 
Note:  1989  has  a  range  from  one  to  three  because  it  is  treated  as  year  one  in  the  study  and 
can  have  no  abnormal  “hits”  for  a  longitudinal  decline. 

Table  4-2 

Summary  of  Abnormal  Readings  for  One  Lung  Function  Test 


Year 

One 

Two 

Three 

Four 

Five 

Six 

1989 

33 

22 

4 

N/A 

N/A 

N/A 

1990 

96 

86 

25 

7 

1 

0 

1991 

68 

46 

15 

6 

2 

0 

1992 

50 

39 

7 

4 

4 

1 

1993 

65 

42 

16 

3 

0 

1 

1994 

65 

29 

7 

2 

1 

0 

1995 

27 

21 

5 

1 

0 

0 
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Related  to  Table  4-2  is  the  information  presented  in  Table  4-3  which  shows  for 
each  lung  function  where  the  abnormal  “hits”  took  place  for  each  of  the  years. 

Table  4-3 

Classification  of  Abnormal  Readings 


Year 

Abnormal 

FVC 

Abnormal 

Ratio 

Abnormal 
Decline  of 
FVC 

Abnormal 
Decline  of 
FEV, 

Abnormal 
Decline  of 
Ratio 

1989 

48 

32 

9 

0 

0 

0 

1990 

146 

116 

38 

27 

24 

25 

1991 

52 

77 

46 

10 

13 

41 

1992 

50 

60 

38 

18 

11 

14 

1993 

54 

63 

48 

18 

8 

24 

1994 

27 

49 

43 

6 

10 

22 

1995 

25 

28 

19 

11 

3 

2 

From  the  information  gathered  by  classifying  each  individual  fimction  test  as 
normal  or  abnormal  based  on  the  six  criteria,  the  next  step  is  creating  control  charts  to 
analyze  which  work  areas  may  be  out-of-control. 


4.3  Control  Charts 

Chart  4-1  is  a  standardized  control  chart  based  on  the  number  of  people  with  an 
abnormal  lung  result  versus  the  total  number  of  people  in  each  of  the  work  areas  for  the 
respective  years.  A  work  area  was  analyzed  only  if  more  than  5  individual  pulmonary 
function  tests  were  administered  for  the  particular  year.  Due  to  the  large  number  of  area- 
year  categories  studied  case  numbers  are  assigned  to  identify  each  area-year  on  the  graph. 
The  interpretation  of  each  case  number  can  be  found  in  appendix  Y.  A  p  value  of  10% 
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was  determined  to  be  used  in  the  application  of  the  control  chart.  This  is  used  because  if 
just  FVC  and  FEVi  were  measured  and  assumed  independent,  9.75%  of  the  normal  tests 
would  be  falsely  classified  as  abnormal. 


Standardized  Control  Chart 
Overall  Population 

;7  =  0.10 


Chart  4-1 

Overall  Population  Standardized  Control  Chart 
All  of  the  points  above  the  3  standard  deviation  Upper  Control  Limit  are 
considered  out-of-control  work-year  areas.  There  are  44  such  areas  of  a  possible  286 
work-year  combinations.  Seven  of  these  areas  include  the  overall  base  average  for  each  of 
the  years  1989-1995.  Table  4-1  above  showed  the  overall  abnormal  percentage  for  each 
of  these  areas;  all  well  above  10%.  The  actual  standard  deviations  above  a  10%  center 
line  ranged  from  3.9  to  13.5  for  the  overall  base  as  displayed  on  the  standardized  control 
chart  above.  A  summary  of  the  44  zone-years  out-of-control  is  shown  in  Table  4-4. 
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Table  4-4 

Out-of-Control  Work  Areas 


Zone 

Years  Out  of  Control 

A1400 

89,  90 

A278A 

90 

A830F 

89,  93 

A867A 

90 

A876C 

90 

Whole  Base 

89,  90,  91,  92,  93,  94,  95 

B36A 

90,  94 

B470B 

91 

B490A 

90,91 

B4B1 

91,92 

B4D 

90 

BSE 

92 

B5J1 

95 

B65A 

90 

B745B 

91 

B770A 

90 

C13D 

93 

C13F 

90 

C13R 

94 

C163A 

90,  95 

C22A 

90 

C22I 

93 

C29A1 

89,  94 

C4020A 

90 

C883A 

94 

C91B1 

94 

K1240 

89,  90,  92,  93,  94 

Due  to  the  fact  the  percentage  of  abnormal  readings  for  nonsmokers  was  much 
lower  than  the  overall  average  for  the  base,  a  standardized  control  chart  for  nonsmokers 
only  also  was  developed.  This  chart  is  shown  in  Chart  4-2.  The  number  of  out-of-control 
work  area-years  in  this  case  dropped  from  44  to  15.  The  overall  base  average  for 
nonsmokers  has  only  one  year  out-of-control  (1990). 
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Chart  4-2 

Nonsmokers  Standardized  Control  Chart 


Table  4-5  shows  the  work  areas  and  given  years  for  which  they  were  out  of 
control  when  smokers  were  excluded  from  the  study. 


Table  4-5 

Nonsmoker  Out-of-Control  Work  Areas 


Zone 

Years  Out  of  Control 

A278A 

90 

A830F 

89, 90,  93 

Whole  Base 

90 

B36A 

90 

B4B1 

91 

B770A 

90 

C13D 

93 

C13F 

90 

C4067A 

93 

K1240 

89, 90,  92,  94 
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Zone  K1240  is  prevalent  in  both  cases.  This  work  zone  is  the  Heat  Plant  located 


in  the  Kittyhawk  area  of  the  base.  A  /?-chart  for  this  specific  zone  is  shown  in  Chart  4-3. 
A/?-chart  of  zone  K1240  also  is  shown  for  the  non-smoking  workers  only  in  Chart  4-4. 


Zone  K1240 

60% 

50% 

40% 

2  30% 

20% 

10% 

0% 

1  2  3  4  5  6  7 

Case  Year 


Chart  4-3 
Zone  K1240 


Zone  K1240 
Nonsmokers 


- ^--P 

^  Pet 

- UCL 

- 1  -Sigma 


Chart  4-4 

Zone  K1240  Nonsmokers 
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As  seen  in  the  above  charts  the  lowest  percentage  of  abnormalities  is  over  20%  (in 
all  seven  years  the  percent  of  abnormalities  are  at  least  one-sigma  above  the  /?-value).  In 
both  cases:  including  and  excluding  the  smoking  population  this  work  area  has  an 
abnormal  rate  which  is  well  above  the  one-sigma  limit  in  all  years  of  the  study.  This  is  a 
zone  which  is  recommended  for  further  study  by  the  74th  SGPMO. 

Montgomery  states  a  process  is  considered  out-of-control  if  2  of  3  consecutive 
points  are  above  the  2-sigma  limits  or  4  of  5  are  above  the  1 -sigma  limit  (20: 117).  Based 
on  this  criteria  other  zones  recommended  for  further  study  based  on  the  overall  population 
are:  B36A  (Heat  Distribution  Center  in  Area  B)  which  had  6  of  the  7  years  at  least  one 
standard  deviation  above  the  /?-value,  B770A  (Heat  Plant  in  Area  B)  which  also  had  6  of 
the  7  years  at  least  one  standard  deviation  above  the  /7-value,  C22I  (88th  Civil  Engineering 
Squadron’s  Project  Painters)  which  had  4  of  the  5  years  at  least  one  standard  deviation 
above  the  /7-value,  and  zone  C163A  (Fire  Stations  #1,  #2  and  #5)  which  had  2  of  the  7 
years  at  least  3  standard  deviations  above  the  /7- value.  Charts  4-5  through  4-12  show  the 
results  of  these  zones  including  and  excluding  the  smoking  population.  Appendix  VI 
reflects  the  Excel  Spreadsheets  used  in  creating  the  charts  for  each  of  the  respective  areas. 
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Chart  4-5 
Zone  B36A 


Zone  B36A 
Nonsmokers 


Chart  4-6 

Zone  B36A  Nonsmokers 


Chart  4-5  shows  zone  B36A  as  having  consistently  high  abnormal  rates.  By 
excluding  people  who  smoke  and  accomplishing  the  p-chart  of  Chart  4-6,  we  see  this 
zone  actually  has  no  abnormalities  for  the  past  2  years  in  which  nonsmokers  were 
measured  (case  year  6  had  zero  nonsmokers  measured). 
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Chart  4-8 

Zone  B770A  Nonsmokers 


Chart  4-7  shows  zone  B770A  as  having  consistently  high  abnormal  rates.  By 
excluding  people  who  smoke  we  see  this  zone  actually  has  no  abnormalities  above  the 
one-sigma  limit  since  case  year  2  of  the  study. 
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Chart  4-9 
Zone  C22I 


Chart  4-9  shows  zone  C22I  as  having  consistently  high  abnormal  rates.  After 
excluding  people  who  smoke,  this  zone  still  has  a  higher  than  1 -sigma  percentage  of 
abnormalities  for  three  of  the  five  years.  Case  year  1  is  1990  and  Case  year  5  is  1995. 
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Chart  4-11 
Zone  Cl 63 A 


Zone  C163A 
Nonsmokers 


1  2  3  4  5  6  7 


—  "  -P 

Pet 

- UCL 

- 1 -Sigma 


Case  Year 


Chart  4-12 

Zone  C163A  Nonsmokers 


Zone  Cl  63  A  is  out  of  control  in  Case  years  2  and  7  with  respect  to  it’s  entire 


population.  However,  when  smokers  are  excluded,  it  is  within  the  three  standard 
deviation  limit  for  Case  year  2  and  the  one  standard  deviation  limit  for  Case  year  7. 


34 


Chapter  5  Summary  and  Recommendations 


5.1  Summary  of  Key  Results 

The  main  result  of  this  study  is  the  74th  SGPMO  now  has  a  tool  to  use  for 
analyzing  the  data  collected  on  the  WPAFB  population.  This  tool  is  the  /?-chart  for  each 
respective  work  area .  To  maintain  this  analysis  they  only  need  to  update  the  results  from 
the  respective  work  areas  each  year  into  a  spreadsheet.  The  initial  results  of  the  study  did 
show  a  large  percentage  of  abnormalities  base  wide  which  indicates  the  population  as  a 
whole  at  WPAFB  is  unhealthy.  Even  with  this  large  percentage  of  abnormalities,  certain 
work  areas  over  the  years  of  the  study  displayed  higher  percentages  of  abnormal  readings 
among  their  workers  than  the  other  work  areas.  These  work  areas  were  K1240,  B36A, 
and  B770A  (all  heat  plants  or  heat  distribution  centers),  C22I  (the  painters  for  the  CE 
squadron),  and  Cl 63 A  (firestations).  The  main  recommendation  from  this  study  is  for  the 
74th  SGPMO  to  recognize  these  particular  work  areas  as  possible  risk  zones. 

Section  5.2  The  Smoker  Issue 

As  stated  above,  the  population  monitored  at  WPAFB  is  on  the  whole  unhealthy. 
Table  5-1  breaks  doAvn  the  total  number  of  individual  measurements  used  in  the  study  by 
smoking  status.  Since  both  of  these  populations  are  binomial  in  nature  (a  test  is  either 
normal  or  abnormal),  in  comparing  the  two  it  is  necessary  to  work  with  the  proportion  of 
abnormalities  and  test  to  see  if  this  is  the  same  (1 1 ;  145).  This  test  is  a  test  of  proportions 
and  the  logic  of  this  test  follows  in  Figure  5-1 . 
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Table  5-1 

Abnormalities  by  Smoking  Status 


Status 

Total  Measurements 

Abnormal  Results 

Percent  Abnormal 

Overall 

4544 

801 

18% 

Nonsmokers 

2090 

287 

14% 

Smokers 

2454 

514 

21% 

•  Null  Hypothesis 

-  Ps  =  Pn 

•  Alternative  Hypothesis 

-  Ps  >  Pn 

•  Test  Statistic  and  Zg  993 

-Z*  =  6.34;  Zo.998=3.09 

•  Conclusion 

-  Reject  the  Null 

Ps  is  the  percent  of  Smokers  with  an  Abnormal  Classification 
Pn  is  the  percent  of  Nonsmokers  with  an  Abnormal  Classification 


Figure  5-1 
Test  of  Proportions 

The  calculations  for  the  above  test  can  be  found  in  appendix  VII.  Since  we  reject 
the  Null  Hypothesis,  we  conclude  there  is  a  true  difference  in  the  two  populations.  It  is 
precisely  for  this  reason  that  the  nonsmoker  population  was  analyzed  on  it’s  own  and 
separate  control  charts  were  constructed  for  capturing  possible  hazardous  work  areas  due 
to  occupational  exposures. 
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5  ■  3  Recommendations  and  Ideas  for  Research 


The  first  recommendation  involves  the  initial  data  entry  of  the  spirometry  results 
into  the  computer.  A  system  should  be  implemented  to  help  reduce  the  number  of  bad 
data  entries.  The  problem  faced  in  this  study  was  the  majority  of  the  ratio  numbers  were 
input  as  FVC  or  FEVi.  There  was  no  way  to  logically  check  the  validity  of  the  FVC  and 
FEVi  entries  once  input  into  the  database  without  going  back  to  the  initial  paper  record 
used  to  input  the  information.  The  current  system  has  the  individual  accomplish  the  exam 
and  a  small  printout  of  the  results  is  generated  by  the  computer;  reference  back  to  Figure 
1-2.  The  information  fi'om  this  printout  is  then  written  on  a  form  which  is  placed  in  the 
person’s  file.  The  information  from  this  form  is  then  read  and  input  to  the  database  by  a 
data  entry  person.  In  each  step  of  the  process  the  potential  for  an  error  exists.  Ideally,  a 
Management  Information  System  should  be  set  up  so  the  results  of  the  exam  are  directly 
placed  into  the  database  fi'om  the  spirometer  itself;  this  will  minimize  human  interaction 
and  thus  decrease  the  error  rates. 

The  second  recommendation  involves  the  classification  of  a  certain  test  as 
abnormal  or  normal.  Instead  of  relying  on  a  population  predicted  value  which  will  classify 
at  least  5%  falsely  as  abnormal,  add  a  box  to  the  spirometry  output  in  which  the  examiner 
can  classify  the  test  as  properly  on  a  “year-to-year”  comparison  basis  or  if  the  individual  is 
tall  and  has  a  genetically  reduced  FEVi  /  FVC  ratio.  This  individual’s  past  test  results 
must  be  accessible  to  the  examiner  for  this  to  work.  This  result  can  then  be  entered  into 
the  database  as  a  binary  variable.  This  will  reduce  a  great  number  of  assumptions  used  in 
this  paper  in  classifying  a  test  as  normal  or  abnormal.  The  main  problem  with  respect  to 
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the  FEVi  /  FVC  ratio  is  there  exists  no  way  for  a  programmer  to  determine  if  it  is  due  to 
the  individual  being  in  great  shape  whereas  the  examiner  would  be  in  a  much  better 
position  to  make  this  determination. 

The  third  recommendation  is  for  a  set  of  regression  coefficients  to  be  developed 
for  the  WPAFB  population  and  to  use  these  as  normal  values.  This  will  help  the  examiner 
if  the  second  recommendation  is  implemented  or  help  the  programmer  feel  more  confident 
about  the  statistical  classification  of  a  test  if  this  recommendation  is  not  implemented. 

The  fourth  recommendation  is  for  the  74th  SGPMO  to  conduct  an  experimental 
design  for  finding  the  true  expected  longitudinal  decline  based  on  age  for  each  of  the  tests. 
Based  on  the  results  of  such  tests,  a  cutoff  for  abnormality  can  be  implemented  to  aid  the 
examiner  in  the  second  recommendation  or  the  programmer  in  the  third  recommendation. 
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Appendix  I:  CONVERT.  SAS  Program 


libname  health  ’user3'; 
run; 

data  health.demo; 
infile  'thesisl.'; 

input  first  $  1  ssan  $  1-9  yr  14-15  sex  $  17  race  $  21  ht  24-25  hyr  29-30  hmo  32-33  hdy  35-36; 
if  index('0 123456789’,first)>0; 
format  htdate  yyinmdd6.; 
htdate  =  mdy(hmo,hdy,hyr); 
drop  hyr  hmo  hdy; 
data  health,  lung; 
infile  'thesis2.'; 

input  first  $  1  ssan  $  1-9  lyr  13-14  Imo  16-17  Idy  19-20  fevl  22-25  fvc  27-30  ratio  32-35  fi^cpred  38-40 
fevpred  45-47  ratpred  52-54; 
ifindexC0123456789',first)>0; 
format  lungdate  yymmdd6.; 
lungdate  =  mdy(lmo,ldy,lyr); 
drop  Imo  Idy  lyr; 
data  health.chem; 
infile  'thesis3.'; 

input  first  $  1  ssan  $  1-9  sgpt  11-13  sgot  16-18  yr  23-24  mo  26-27  dy  29-30  ap  32-33  ggt  37-39 
bili  42-45  albumin  52-55; 
if  index(‘0 123456789\first)>0; 
format  chemdate  yynundd6.; 
chemdate  =  mdy(mo,dy,yr); 
drop  mo  dy  yr; 
data  health.blood; 
infile  'thesis4.'; 

input  first  $  1  ssan  $  1-9  wbc  $  11-14  yr  20-21  mo  23-24  dy  26-27  hemcrit  29-33; 
ifindexC0123456789’,first)>0; 
format  blddate  yymmdd6.; 
blddate  =  mdy(mo,dy,yr); 
drop  mo  dy  yr; 
data  health.zone; 
infile  'thesis5.'; 

input  first  $  1  ssan  $  1-9  zone  $  11-18  syr  22-23  smo  25-26sdy  28-29  eyr  33-34  emo  36-37  edy  39-40; 
ifindexC0123456789',first)>0; 
format  stdate  yymmdd6.  enddate  yymmdd6.; 
stdate  =  mdy(smo,sdy,syr); 
enddate  =  mdy(emo,edy,eyr); 
drop  smo  sdy  syr  emo  edy  eyr; 
data  health.history; 
infile  'thesis6,'; 

input  first  $  1  ssan  $  1-9  blddis  $  11  asthma  $  15  lung  $  19  liver  $  23  heptitis  $  27  jaundice  $  31  smkl  35 
smk2  39-40  smk3  43-44; 
ifindexC0123456789',first)>0; 
data  health.drinking; 
infile  'thesis7.'; 

input  first  $  1  ssan  $  1-9  liquor  11-12  beer  15-16  wine  19-20; 
if  index(*0 123456789’, first)>0; 
run; 


39 


Appendix  II:  FVCRAW.SAS  Program 


libname  health  'user3*; 
run; 

options  Is  =  75  ; 
proc  sort  data  =  health,  lung; 
by  ssan; 
run; 

data  _null_; 
set  health.lung; 
by  ssan; 

file  print  notitles; 
if  first.ssan  then  do; 
put@l  ssan@llfvc  @; 
n  =  16; 
end; 

if  first.ssan  =  0  and  last.ssan  =  0  then  do; 
put@nfvc  @; 
n  =  n+5; 
retain  n; 
end; 

if  last  ssan  then  do; 
put@nfvc  @75  first; 
end; 
run; 
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Appendix  HI:  MERGEALL.SAS  Program 


data  a; 

infile  'demo.raw'; 

input  ssan  $  1-9  dobyr  11-12  sex  $  14  race  $  15 
ht  16-17; 
run; 

proc  sort  data  =  a;  by  ssan ;  run; 
datab; 

infile  'fevl.raw'; 

input  ssan  $  1-9  fevll  11-15  fevl2  16-20  fevl3 
21-25  fevl4  26-30 

fevl5  31-35  fevl6  36-40  fevl7  41-45  fevl8 
46-50; 
run; 

proc  sort  data  =  b;  by  ssan ;  run; 
datac; 

infile  'fi^c.raw'; 

input  ssan  $  1-9  fvcl  1 1-15  fvc2  16-20  fvc3 
21-25  fvc4  26-30 

fvc5  31-35  fvc6  36-40  fvc7  41-45  lvc8 
46-50; 
run; 

proc  sort  data  =  c;  by  ssan  ;  run; 
datad; 

infile  'ratio.raw'; 

input  ssan  $  1-9  ratiol  11-15  ratio2  16-20 
ratio3  21-25  ratio4  26-30 

ratio5  31-35  ratio6  36-40  ratio7  41-45  ratio8 
46-50; 
run; 

proc  sort  data  =  d;  by  ssan  ;  run; 
datae; 

infile  'fi^cpred.raw'; 

input  ssan  $  1-9  fvcpdl  11-15  fvcpd2  16-20 
fvcpd3  21-25  lVcpd4  26-30 

lVcpd5  31-35  fvcpd6  36-40  fVcpd7  41-45 
fvcpdS  46-50; 
run; 

proc  sort  data  =  e;  by  ssan  ;  run; 
dataf; 

infile  'fevpred.raw'; 

input  ssan  $  1-9  fevpdl  11-15  fevpd2  16-20 
fevpd3  21-25  fevpd4  26-30 

fevpd5  31-35  fevpd6  36-40  fevpd7  41-45 
fevpd8  46-50; 


run; 

proc  sort  data  =  f;  by  ssan ;  run; 
datag; 

infile  'ratpred.raw'; 

input  ssan  $  1-9  ratpdl  11-15  ratpd2  16-20 
ratpd3  21-25  ratpd4  26-30 

ratpd5  31-35  ratpd6  36-40  ratpd7  41-45 
ratpdS  46-50; 
run; 

proc  sort  data  =  g;  by  ssan ;  run; 
datah; 

infile  'lungdate.raw'; 

input  ssan  $  1-9  yrl  11-12  mol  13-14  dyl  15- 
16  yr2  18-19  mo2  20-21 

dy2  22-23  yr3  25-26  mo3  27-28  dy3  29-30 
yr4  32-33  mo4  34-35 

dy4  36-37  yr5  39-40  mo5  41-42  dy5  43-44 
yr6  46-47  mo6  48-49 

dy6  50-51  ytl  53-54  mo7  54-56  dy7  57-58 
yr8  60-61  mo8  62-63 
dy8  64-65; 

format  Idtl  ldt2  ldt3  ldt4  ldt5  ldt6  ldt7  IdtS 
yymmdd6.; 

Idtl  =  mdy(mol,dyl,yrl); 
ldt2  =  mdy(mo2,dy2,yr2); 
ldt3  =  mdy(mo3,dy3,yr3); 
ldt4  =  mdy(mo4,dy4,yr4); 
ldt5  =  mdy(mo5,dy5,yr5); 
ldt6  =  mdy(mo6,dy6,yr6); 
ldt7  =  mdy(mo7,dy7,yr7); 
ldt8  =  mdy(mo8,dy8,yr8); 
drop  yrl  mol  dyl  yr2  mo2  dy2  yr3  mo3  dy3  yr4 
mo4  dy4  yr5  mo5  dy5  yr6 

mo6  dy6  yrl  mo7  dy7  yr8  mo8  dy8; 
run; 

proc  sort  data  =  h;  by  ssan  ;  run; 
data  i; 

infile  'sgpt.raw'; 

input  ssan  $  1-9  sgptl  11-14  sgpt2  15-18 
sgpt3  19-22  sgpt4  23-26 

sgpt5  27-30  sgpt6  31-34  sgpt7  35-38  sgpt8 
39-42; 
run; 

proc  sort  data  =  i;  by  ssan  ;  nm; 
dataj; 
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infile  'sgot.raw'; 

input  ssan  $  1-9  sgotl  11-14  sgot2  15-18  sgotS 
19-22  sgot4  23-26 

sgot5  27-30  sgot6  31-34  sgot7  35-38  sgotS 
39-42; 
run; 

proc  sort  data  =  j;  by  ssan  ;  run; 
datak; 

infile  'ap.raw'  ; 

input  ssan  $  1-9  apl  1 1-14  ap2  15-18  ap3 
19-22  ap4  23-26 

ap5  27-30  ap6  31-34  ap7  35-38  ap8 
39-42; 
run; 

proc  sort  data  =  k;  by  ssan ;  run; 
datal; 

infile  'ggt.raw' ; 

input  ssan  $  1-9  ggtl  11-14  ggt2  15-18  ggt3 
19-22  ggt4  23-26 

ggt5  27-30  ggt6  31-34ggt7  35-38  ggt8 
39-42; 
run; 

proc  sort  data  =  1;  by  ssan ;  run; 
datam; 

infile  "bili.raw'; 

input  ssan  $  1-9  bilil  11-15  bili2  16-20  bili3 
21-25  bili4  26-30 

bili5  31-35  bili6  36-40  bili7  41-45  bili8 
46-51; 
run; 

proc  sort  data  =  m;  by  ssan ;  run; 
data  n; 

infile  'albumin.raw' ; 

input  ssan  $  1-9  albuminl  11-15  albunun2  16- 
20  albuniin3  21-25 

albumin4  26-30  albuinin5  31-35  albumin6 
36-40  albumin7  41-45 
albunun8  46-50; 
run; 

proc  sort  data  =  n;  by  ssan ;  run; 
data  o; 

infile  'chemdate.raw'; 

input  ssan  $  1-9  yrl  11-12  mol  13-14  dyl  15- 
16  yr2  18-19  mo2  20-21 

dy2  22-23  yr3  26-27  mo3  28-29  dy3  30-31 
yr4  34-35  mo4  36-37 


<fy4  38-39  yr5  42-43  mo5  44-45  dy5  46-47 
yr6  50-51  mo6  52-53 

dy6  54-55  yr7  58-59  mo7  60-61  dy7  62-63 
yr8  66-67  mo8  68-69 
dy8  70-71; 

format  cdtl  cdt2  cdt3  cdt4  cdt5  cdt6  cdt7  cdt8 
yymmdd6.; 

cdtl  =  mdy(mol,dyl,yrl); 
cdt2  =  mdy(mo2,^2,yr2); 
cdt3  =  mdy(mo3,dy3,yr3); 
cdt4  =  mdy(mo4,dy4,yr4); 
cdt5  =  mdy(mo5,dy5,yr5); 
cdt6  =  mdy(mo6,dy6,yr6); 
cdt7  ==  mdy(mo7,dy7,yr7); 
cdt8  =  mdy(mo8,dy8,yr8); 
drop  yrl  mol  dyl  yr2  mo2  dy2  yr3  mo3  dy3  yr4 
mo4  dy4  yr5  mo5  (fy5  yr6 

mo6  dy6  yr7  mo7  dy7  yr8  mo8  dy8; 
run; 

proc  sort  data  =  o;  by  ssan  ;  run; 
datap; 

infile  "blddate-raw' ; 

input  ssan  $  1-9  yrl  11-12  mol  13-14  dyl  15- 
16  yr2  18-19  mo2  20-21 

dy2  22-23  yr3  26-27  mo3  28-29  dy3  30-31 
yr4  34-35  mo4  36-37 

dy4  38-39  yr5  42-43  mo5  44-45  dy5  46-47 
yr6  50-51  mo6  52-53 

dy6  54-55  yr7  58-59  mo7  60-61  dy7  62-63 
yr8  66-67  mo8  68-69 
dy8  70-71; 

format  bdtl  bdt2  bdt3  bdt4  bdt5  bdt6  bdt7  bdt8 
yymmdd6.; 

bdtl  =mdy(mol,dyl,yrl); 

bdt2  =  mdy(mo2,dy2,yr2); 

bdt3  =  mdy(mo3,dy3,yr3); 

bdt4  =  mdy(mo4,dy4,yr4); 

bdt5  =  mdy(mo5,dy5,yr5); 

bdt6  =  mdy(mo6,dy6,yr6); 

bdt7  =  mdy(mo7,dy7,yr7); 

bdt8  =  mdy(mo8,dy8,yr8); 

drop  yrl  mol  dyl  yr2  mo2  dy2  yr3  mo3  dy3  yr4 

mo4  dy4  yr5  mo5  dy5  yr6 

mo6  dy6  yr7  mo7  dy7  yr8  mo8  dy8; 
run; 

proc  sort  data  =  p;  by  ssan  ;  run; 
data  q; 

infile  'wbc.raw' ; 

input  ssan  $  1-9  wbcl  11-17  wbc2  18-25 
wbc3  26-33  wbc4  34-41 
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wbc5  42-49  wbc6  50-57  wbc7  58-65 
wbc8  66-74; 
run; 

proc  sort  data  =  q;  by  ssan ;  run; 
datar; 

infile  'hemcrit.raw'; 

input  ssan  $  1-9  heml  11-17  hem2  18-25 
hem3  26-33  hem4  34-41 

hem5  42-49  hem6  50-57  hem7  58-65 
hem8  66-74; 
run; 

proc  sort  data  =  r;  by  ssan  ;  run; 
datas; 

infile  'zone.raw'  ; 

input  ssan  $  1-9  zonel  $11-19  zone2  $20-28 
zone3  $29-37  zone4  $38-46 

zone5  $47-55  zone6  $56-64  zone7  $65-73  ; 

run; 

proc  sort  data  =  s;  by  ssan  ;  run; 
datat; 

infile  'stdate.raw'  ; 

input  ssan  $  1-9  yrl  11-12  mol  13-14  dyl  15- 
16  yr2  20-21  mo2  22-23 

dy2  24-25  yr3  29-30  mo3  31-32  dy3  33-34 
yr4  38-39  mo4  40-41 

dy4  42-43  yr5  47-48  mo5  49-50  dy5  51-52 
yr6  56-57  mo6  58-59 

dy6  60-61  yrl  65-66  mo7  67-68  dy7  69-70; 
format  sdtl  sdt2  sdt3  sdt4  sdt5  sdt6  sdt7 
yymmdd6.; 

sdtl  =  mdy(mol,dyl,yrl); 
sdt2  =  mdy(mo2,dy2,yT2); 
sdt3  =  mdy(mo3,dy3,yr3); 
sdt4  =  mdy(mo4,dy4,yr4); 
sdt5  =  mdy(mo5,dy5,yr5); 
sdt6  =  mdy(mo6,dy6,yr6); 
sdt7  =  mdy(mo7,dy7,yr7); 
drop  yrl  mol  dyl  yr2  mo2  dy2  yr3  mo3  dy3  yr4 
mo4  dy4  yr5  mo5  dy5  yr6 
mo6  dy6  yrl  mo7  dy7; 
run; 

proc  sort  data  =  t;  by  ssan  ;  run; 
datau; 

infile  'enddate.raw' ; 

input  ssan  $  1-9  yrl  11-12  mol  13-14  dyl  15- 
16  yr2  20-21  mo2  22-23 

dy2  24-25  yr3  29-30  mo3  31-32  dy3  33-34 
yr4  38-39  mo4  40-41 


dy4  42-43  yr5  47-48  mo5  49-50  dy5  51-52 
yr6  56-57  mo6  58-59 

dy6  60-61  yrl  65-66  mo7  67-68  dy7  69-70; 
format  edtl  edt2  edt3  edt4  edt5  edt6  edt7 
yymmdd6.; 

edtl  =  mdy(mol,dyl,yrl); 

edt2  =  mdy(mo2,dy2,yr2); 
edt3  =  mdy(mo3,dy3,yr3); 
edt4  =  mdy(mo4,dy4,yr4); 
edt5  =  mdy(mo5,dy5,yr5); 
edt6  =  mdy(mo6,dy6,yr6); 
edt7  =  mdy(mo7,dy7,yr7); 
drop  yrl  mol  dyl  yr2  mo2  dy2  yr3  mo3  dy3  yr4 
mo4  dy4  yr5  mo5  dy5  yr6 
mo6  dy6  yrl  mo7  dy7; 
run; 

proc  sort  data  =  u;  by  ssan  ;  run; 
data  v; 

infile  'history. raw’; 

input  ssan  $  1-9  blddis  $  11  asthma  $  12  lung  $ 
13  liver  $  14 

heptitis  $  15  jaundice  $  16  smkl  17-18 
smk2  19-21  smk3  22-23; 
run; 

proc  sort  data  =  v;  by  ssan  ;  run; 
data  w; 

infile  'drinking.raw'; 

input  ssan  $  1-9  liquor  11-13  beer  14-16  wine 

17-19; 

run; 

proc  sort  data  =  w;  by  ssan  ;  run; 


libname  health  'user3'; 
data  health.wpafbl; 

merge  a  b  (in=inl)  c  d  e  f  g  h  i  (in=in2)  J  k  1  m 
n  0  p  q  (in=in3) 
r  s  t  u  V  w; 
by  ssan; 

if  ini  or  in2  or  in3; 
if  zonel  =  '  '  then  delete; 
run; 

proc  contents; 
run; 
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Appendix  IV:  LUNGYRSl  SAS  and  LUNGYRS2.SAS  Programs 


Section  IV.  1  LUNGYRS 1 .  SAS  Program 


libname  health  'user3'; 
run; 

data  a  yl987  yl988  yl989  yl990  yl991  yl992 
yl993  yl994  yl995; 

set  health. wpafb2  (keep  =  ldtl-ldt8  fvcl-fvc8 
fevll-fevl8  ssan 

ratiol-ratio8  fVq)dl-fvq)d8  fevpdl-fevpdS 
ratpdl-ratpd8 

zone  1 -zone?  sdtl-sdt?  edtl-edt?); 
array  fvc[*]  fvcl-fvc8; 
array  fevl[*]  fevl  l-fevl8; 
array  ratio[*]  ratiol-ratioS; 
array  ratpd[*]  ratpdl-ratpd8; 
array  fvq)d[*]  fvq)dl-^cpd8; 
array  fevpd[*]  fevpdl-fevpd8; 
array  ldt[*]  ldtl-ldt8; 
array  ldtyr[*]  ldtyrl-ldtyr8; 
array  keep[*]  keepl-keep8; 
array  zone[*]  zonel-zoneS; 
array  sdt[*]  sdtl-sdt8; 
array  edt[*]  edtl-edt8; 
zone8  = ' 
sdt8  = 
edt8  = 

format  cratio  6.2; 

/* 

Create  Proper  Work  Areas 

*! 

do  j  =  1  to  7; 

if  zone[j]  =  'A1400'  then  zone[j]  = 
substr(zone|j],l,6); 

else  if  zone[j]  =  'B6000'  then  zone|j]  = 
substr(zone[j],l,6); 

else  if  zoneU]  =  'C4020'  then  zoneQ]  = 
substr(zone[j],l,6); 

else  if  zone[j]  =  'K1084'  then  zoneU]  = 
substr(zone[j],l,6); 
else  if 

index('B70.B18.B22.B36.B63.B79.C89’,substr( 
zone[j],l,3))>  0 

then  zonelJ]  =  substr(zone[j],l,4); 
else  if 

index('B40.C13.C19.C22.C28.C59.C70’,substr( 

zone[j],l,3))>0 

then  zone|j]  =  substr(zone|j],l,4); 
else  if  substr(zone|j],l,4)  =  'B24A'  then 
zoneU]  =  'B24A'; 

else  if  substr(zoneU],l,4)  =  'B24B'  then 
zoneU]  =  'B24B'; 


else  if  substr(zoneU],l,4)  =  'B45B'  then 
zoneU]  =  WSB’; 

else  if  substr(zoneU],l»4)  =  "BbSA'  then 
zoneU]  =  'B65A'; 

else  if  substr(zoneU],l,3)  =  ’B4D'  then  zoneU] 
=  'B4D'; 

else  if  substr(zoneU],l>3)  =  'B4E'  then  zoneU] 
=  'B4E’; 

else  if  substr(zoneU],l,3)  =  'BSC  then  zoneU] 
=  'B5C'; 

else  if  substr(zoneU],l73)  =  "BSO'  then  zoneU] 
=  'B5D'; 

else  if  substr(zoneU],l,3)  =  'B5E'  then  zoneU] 
=  'B5E'; 

else  if  substr(zoneU],E3)  =  "BSH'  then  zoneU] 
=  'B5H'; 

else  if  substr(zoneU],l,2)  =  'CE'  then  zoneU] 

=  'CE'; 

else  if  substr(zoneU],l,3)  =  'C40' 
then  zoneU]  =  substr(zoneU],l,6); 
else  zoneU]  =  substr(zoneU],l,5); 
end; 

doj  =  1  to7; 

if  sdtU]  .  and  edtU]  =  .  then  edtU]  = 
’31dec99'd; 
end; 

doj  =  1  to  8; 

if  fvcU]  .  then  keepU]  =  1; 
end; 

if  keepl  =  1  or  keep2  =  1  or  keep3  =  1  or  keep4 
=  1  or  keeps  =  1 

or  keep6  =  1  or  keep?  =  1  or  keep8  =  1; 
doj  =  1  to  8; 

if  keepU]  =  .  then  keepU]  =  0; 
if  keepl  =  keep2  =  1  and  Idtl  =  .  then  Idtl  = 
ldt2-365; 
end; 

doj  =  1  to  7; 

if  keepU]  =  keepU+1]  =  1  then  do; 
if  IdtU]  =  ldtU+1]  then  do; 
ldtU+1]  = .; 
keepU+1]  =  0; 
end; 
end; 
end; 

doj  =  1  to  8; 

if  keepU]  =  1  then  IdtyrU]  =  year(ldtU]); 
if  IdtyrU]  .  then  do; 
if  IdtyrU]  <  1988  then  do; 
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vc=  fvclj]; 
fev  =  fevl|j]; 
fev_fvc  =  ratioU]; 
predfvc  =  fvqpd|j]; 
predfev  =  fevpdjj]; 
prdratio=  ratpd[j]; 
tyear  =  IdtyrU]; 

if  VC  .  and  fev  .  then  cratio 

=(fev/vc)*100; 

if  sdt[j]  <=  ldt[j]  <=  edt[j]  then  czone  = 
zone[j]; 

if  substr(czone.l,l)  =  ' '  then  czone  =  zonel; 

output  yl987; 
end; 

else  if  Idtyrjj]  ==  1988  then  do; 
vc=  fv^clj]; 
fev  =  fevlO]; 
fevfvc  =  ratioQ]; 
predfvc  =  fv'cpdlj]; 
predfev  =  fevpd[j]; 
prdratio=  ratpdUJ; 
tyear  =  ldtyr[j]; 

if  VC  .  and  fev  .  then  cratio 
=(fev/vc)*100; 

if  sdt[j]  <=  ldt[j]  <=  edt[j]  then  czone  = 
zone[j]; 

if  substr(czone,l,l)  =  ’  ’  then  czone  =  zonel; 

output  yl988; 
end; 

else  if  ldtyr[j]  =  1989  then  do; 
vc=  fv^cU]; 
fev  =  fevl[j]; 
fev_fvc  =  ratio[j]; 
predfvc  =  fv^cpdjj]; 
predfev  =  fevpdjj]; 
prdratio=  ratpdjjj; 
tyear  =  ldtyr[j]; 

if  VC  .  and  fev  ,  then  cratio 
=(fev/vc)*100; 

if  sdtjj]  <=  ldt[j]  <=  edtjj]  then  czone  = 
zone[j]; 

if  substr(czone,l,l)  =  ' '  then  czone  =  zonel; 

output  yl989; 
end; 

else  if  Idtyrjj]  ==  1990  then  do; 

VC  =  fvclj]; 
fev  =  fevl[j]; 
fev_fvc  =  ratiolJ]; 
predfvc  =  fv^cpd[j]; 
predfev  =  fevpdjj]; 
prdratio=  ratpdjj]; 
tyear  =  Idtyrjj]; 


if  vc  ,  and  fev  .  then  cratio 
=(fev/vc)*100; 

if  sdtjj]  <=  Idtjj]  <=  edt[j]  then  czone  = 
zonejj]; 

if  substr(czone,l,l)  = ' '  then  czone  =  zonel; 

output  y 1990; 
end; 

else  if  Idtyrjj]  =  1991  then  do; 
vc=  fv^cjj]; 
fev  =  fevljj]; 
fev_fvc  =  ratiojj]; 
predfVc  =  fvcpdjj]; 
predfev  =  fevpdjj]; 
prdratio=  rat^jj]; 
tyear  =  Idtyrjj]; 

if  VC  .  and  fev  ^=  .  then  cratio 
=(fev/vc)*100; 

if  sdt[j]  <=  Idtjj]  <=  edtO]  then  czone  = 
zonejj]; 

if  substr(czone,l,l)  = ' '  then  czone  =  zonel; 

output  yl991; 
end; 

else  if  Idtyrjj]  =  1992  then  do; 
vc=  fvcjj]; 
fev  =  fevljj]; 
fev_fvc  =  ratiojj]; 
predfVc  =  fvcpdjj]; 
predfev  =  fevpdjj]; 
prdratio=  ratpdjj]; 
tyear  =  Idtyrjj]; 

if  VC  .  and  fev  .  then  cratio 
=(fev/vc)*100; 

if  sdtjj]  <=  Idtjj]  <=  edtjj]  then  czone  = 
zonejj]; 

if  substr(czone,l,l)  =  ' '  then  czone  =  zonel; 

output  yl992; 
end; 

else  if  Idtyrjj]  =  1993  then  do; 
vc=  fvcjj]; 
fev  =  fevl|j]; 
fev_fvc  =  ratiojj]; 
predfvc  =  fvcpdjj]; 
predfev  =  fevpdjj]; 
prdratio=  ratpdjj]; 
tyear  =  Idtyrjj]; 

if  VC  .  and  fev  .  then  cratio 
=(fev/vc)*100; 

if  sdtjj]  <=  Idtjj]  <=  edtjj]  then  czone  = 
zonejj]; 

if  substr(czone,l,l)  =  * '  then  czone  =  zonel; 

output  yl993; 
end; 
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else  if  ldtyr[j]  =  1994  then  do; 
vc=  fvc\j]; 
fev  =  fevl[j]; 
fev_fvc  =  ratiolj]; 
predfvc  =  fv^qjdjj]; 
predfev  =  fevpdU]; 
prdratio=  ratpdfj]; 
tyear  =  Idtyrlj]; 

if  VC  .  and  fev  .  then  cratio 
=(fev/vc)*100; 

if  sdt[j]  <=  ldt[j]  <=  edt[j]  then  czone  = 
zone[j]; 

if  substr(czone,l,l)  =  ’  ’  then  czone  =  zonel; 

output  yl994; 
end; 

else  if  Idtyrjj]  =  1995  then  do; 
vc=  fv^clj]; 
fev  =  fevlU]; 
fev_fvc  =  ratioOl; 
predfvc  =  fv^cpdjj]; 
predfev  =  fevpd|j]; 
prdratio=  ratpdjj]; 
tyear  =  ldtyr[j]; 

if  VC  .  and  fev  .  then  cratio 
=(fev/vc)*100; 

if  sdtjj]  <=  IdtO]  <=  edt|j]  then  czone  = 
zone[j]; 

if  substr(czone,l,l)  =  '  ’  then  czone  =  zonel; 

output  y 1995; 
end; 
end; 
end; 
run; 

datayl987; 

set  yl987  (keep  =  ssan  vc  fev  fev  fVc 
predfvc  predfev  prdratio 
cratio  czone); 

vc87  =  vc;  fev87  =  fev;  ratio87  =  fev_fvc; 
pfvc87  =  predfvc; 

pfev87  =  predfev;  pratio87  =  prdratio;  cratio87 
=  cratio; 

czone87  =  czone; 
id  =  ssan|r_87'; 
run; 

datayl988; 

set  yl988  (keep  =  ssan  vc  fev  fevjfVc  predfvc 
predfev  prdratio 

cratio  czone); 

vc88  =  vc;  fev88  =  fev;  ratio88  =  fev_fvc; 
pfvc88  =  predfVc; 

pfev88  =  predfev;  pratio88  =  prdratio;  cratio88 
=  cratio; 


czone88  =  czone; 
id  =  ssan||'_88'; 
run; 

data  yl989; 

set  yl989  (keep  =  ssan  vc  fev  fev  fvc 
predfvc  predfev  prdratio 
cratio  czone); 

vc89  =  vc;  fev89  =  fev;  ratio89  =  fev__fvc; 
pfvc89  =  predfvc; 

pfev89  =  predfev;  pratio89  =  prdratio;  cratio89 
=  cratio; 

czone89  =  czone; 
id  =  ssan||'_89’; 
run; 

datayl990; 

set  yl990  (keep  =  ssan  vc  fev  fev  fvc  predfvc 
predfev  prdratio 

cratio  czone); 

vc90  =  vc;  fev90  =  fev;  ratio90  =  fev  fvc; 
pfvc90  =  predfvc; 

pfev90  =  predfev;  pratio90  =  prdratio;  cratio90 
=  cratio; 

czone90  =  czone; 
id  =  ssan||'_90’; 
run; 

datayl991; 

set  y  1 99 1  (keep  =  ssan  vc  fev  fev  fvc 
predfVc  predfev  prdratio 
cratio  czone); 

vc91  =  vc;  fev91  =  fev;  ratio91  ==  fev  fvc; 
pfvc91  =  predfvc; 

pfev91  =  predfev;  pratio91  =  prdratio;  cratio91 
=  cratio; 

czone91  =  czone; 
id  =  ssan|i’_9r; 
run; 

datayi992; 

set  yl992  (keep  =  ssan  vc  fev  fev_fVc  predfvc 
predfev  prdratio 

cratio  czone); 

vc92  =  vc;  fev92  =  fev;  ratio92  =  fev  fvc; 
pfvc92  =  predfvc; 

pfev92  =  predfev;  pratio92  ==  prdratio;  cratio92 
=  cratio; 

czone92  ==  czone; 
id  =  ssan||'_92'; 
run; 

datayl993; 

set  yl993  (keep  =  ssan  vc  fev  fev_fvc 
predfvc  predfev  prdratio 
cratio  czone); 
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vc93  -  vc;  fev93  =  fev;  ratio93  =  fev  fvc; 
pfvc93  =  predfvc; 

pfev93  =  predfev;  pratio93  =  prdratio;  cratio93 
=  cratio; 

czone93  =  czone; 
id  =  ssan||'_93’; 
run; 

datayl994; 

set  yl994  (keep  =  ssan  vc  fev  fev  fvc  predfvc 
predfev  prdratio 

cratio  czone); 

vc94  =  vc;  fev94  =  fev;  ratio94  =  fev_fvc; 
pfvc94  =  predfvc; 

pfev94  =  predfev;  pratio94  =  prdratio;  cratio94 
=  cratio; 

czone94  =  czone; 
id  =  ssan||'_94'; 
run; 

datayl995; 


set  y 1 995  (keep  =  ssan  vc  fev  fev_fvc 
predfvc  predfev  prdratio 
cratio  czone); 

vc95  =  vc;  fev95  =  fev;  ratio95  =  fev_fvc; 
pfvc95  =  predfvc; 

pfev95  =  predfev;  pratio95  =  prdratio;  cratio95 
=  cratio; 

czone95  =  czone; 
id  =  ssan||'_95'; 
run; 

data  health,  lung  (drop  =  vc  fev  fev  fvc  predfvc 
predfev  prdratio  cratio 
czone); 

merge  yl987  yl988  yl989  yl990  yl991  yl992 
yl993  yl994  yl995; 
by  ssan; 
run; 

proc  contents  data  =  health.lung; 
run; 


Section  IV.2  LUNGYRS2.SAS  Proeram 


libname  health  'user3‘; 

run; 

data  a; 

set  health.  Iung3; 
array  cratio[*]  cratio87-cratio95; 
array  czone[*]  czone87-czone95; 
array  fev[*]  fev87-fev95; 
array  pfev[*]  pfev87-pfev95; 
array  pfvc[*]  pfvc87-pfvc95; 
array  pratio[*]  pratio87-pratio95; 
array  ratio[*]  ratio87-ratio95; 
array  abfvc[*]  abfVcl-abfvc9; 
array  abfev[*]  abfevl-abfev9; 
array  abrat[*]  abratl-abrat9; 
array  ablfvc[*]  ablfVcl-ablfvc9; 
array  ablfev[*]  ablfevl-ablfev9; 
array  ablrat[*]  ablratl-ablrat9; 
array  vc[*]  vc87-vc95; 
array  keep[*]  keepl-keep9; 
array  pval[*]  pvall-pval9; 
array  pvala[*]  pvalal"pvala9; 
array  smoke[*]  smokel-smoke9; 
do  j  =  1  to  9  ; 

if  vc|j]  .  then  keep|j]  =  1;  else  keepO]  =  0; 

if  keep[j]  =  1  then  smoke[j]  =  smoker;  else 
smoke[j]  =  0; 
end; 

/* 

Define  Abnormalities  based  on  Percent 
Predicted  Values  of  FVC,  FEVl 


and  the  Observed  Value  of  the  Ratio 
FEVl/FVC 
*/ 

if  dobyr  >  .  and  ht  >  .  then  do; 
htcm  =  2.54*ht; 
do  j  =  1  to  9; 

age  =  (j+86)-dobyr; 
if  sex  =  'M'  then  do; 
if  age  <  25  then  do; 

pval[j]  =  -4.808  +  .045*age  +  .046*htcm; 
end; 
else  do; 

pval|j]  =  -4.203  -  .027*age  +  .052*htcm; 
end; 
end; 

else  if  sex  =  'F*  then  do; 
if  age  <  25  then  do; 

pvallj]  =  -2.703  +  .085*age  +  .027*htcm; 

end; 

else  do; 

pval[j]  =  -0.794  -  .021*age  +  .027*htcm; 

end; 

end; 

pfevLj]  =  (fev[j]/pvallj])*100; 
end; 

doj  =  1  to  9; 

age  =  (j+86)-dobyr; 
if  sex  =  'M'  then  do; 
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if  age  <  25  then  do; 

pvalaO]  =  -5.508  +  .078*age  +  .05*htcm; 

end; 

else  do; 

pvala[j]  =  -5.459  -  .029*age  +  .065*htcm; 
end; 
end; 

else  if  sex  = 'F'  then  do; 
if  age  <  20  then  do; 

pvalajj]  =  -3.469  +  .092*age  +  .033*htcni; 
end; 
else  do; 

pvala[j]  =  -1.774  -  .022*age  +  .037*htcm; 
end; 
end; 

pfvc[j]  =  (vc[j]/pvala[j])*100; 
end; 
end; 

if  dobyr  =  .  then  dobyr  =  50; 
doj  =  1  to  9; 
if  keeplj]  =  1  then  do; 
if  sex  =  'M'  or  sex  =  ' '  then  do; 
if  ((86+j)  -  dobyr)  >39  then  do; 
if  pfvc|j]  .  then  do; 

if  p^c[j]  <  73.4  then  abfvclJ]  =  1;  else 
abfvc[j]  =  0; 
end; 

else  abfvc|j]  =  0; 
if  pfev[j]  .  then  do; 

if  pfev|j]  <  77.2  then  abfevjj]  =  1;  else 
abfev[j]  =  0; 
end; 

else  abfevlJ]  =  0; 
end; 
else  do; 

if  pfvc[j]  '^= .  then  do; 
if  pfvc[j]  <81.1  then  abfvc[j]  =  1;  else 
abfvcIJ]  =  0; 
end; 

else  abfvclj]  =  0; 
if  pfev[j]  .  then  do; 
if  pfev|j]  <  79.1  then  abfev|j]  =  1;  else 
abfev[j]  =  0; 
end; 

else  abfevlJ]  =  0; 
end; 
end; 

else  if  sex  =  'F'  then  do; 
if  ((86+j)  -  dobyr)  >39  then  do; 
if  pfvcjj]  *= .  then  do; 


if  pfvcO]  <  75.2  then  abfvc|j]  =  1;  else 
abfvclj]  =  0; 
end; 

else  abfvclj]  =  0; 
if  pfevlj]  .  then  do; 
if  pfev[j]  <  77.9  then  abfev|j]  =  1;  else 
abfevlJ]  =  0; 
end; 

else  abfevlJ]  =  0; 
end; 
else  do; 

if  pfvcU]  ''= .  then  do; 
if  pfvc|j]  <  76.9  then  abfvc[j]  =  1;  else 
abfvclj]  =  0; 
end; 

else  abfvclj]  =  0; 
if  pfev[j]  ''= .  then  do; 
if  pfev[j]  <  70.3  then  abfev[j]  =  1;  else 
abfevlj]  =  0; 
end; 

else  abfevlj]  =  0; 
end; 
end; 

if  cratiolj]  .  then  do; 

if  cratiolj]  <  70  then  abratlj]  =  1;  else  abratlj] 
=  0; 
end; 

else  abratlj]  =  0; 
end; 
end; 

/* 

Define  abnormalities  based  on  longitudinal 
data.  Based  on 

Percentage  decrease  of  observed  FVC,  FEVl, 
and  Ratio 
*/ 

doj  =  1  to  8; 

if  keeplj]  =  keeplj+1]  =  1  then  do; 
if  vc^]  >  0  and  vclj+1]  >  0  and  vclj]  >  vclj+1] 
then  do; 

if  ((vclj]-vclj+l])/vclj])  >  .15  then 
ablfvclj+1]  =  1; 

else  ablh'clj+l]  =  0; 
end; 

if  fevlj]  >  0  and  fevg+l]  >  0  and  fevlj]  > 
fevlj+1]  then  do; 

if  ((fevlj]-fevlj+l])/fevlj])  >  .15  then 
ablfevlj+l]  =  1; 

else  ablfevO+l]  =  0; 
end; 

if  cratiolj]  >  0  and  cratiolj+1]  >  0  and  cratiolj] 
>  cratiolj+1] 
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then  do; 

if  ((cratio[j]-cratio[j+l])/cratio|j])  >  .10  then 
ablrat|J+l]  =  1; 

else  ablrat[j+ll  =  0; 
end; 
end; 
else  do; 

ablfvcD+1]  =  0; 
ablfev[j+l]  =  0; 
ablrat[j+l]  =  0; 
end; 
end; 

doj  =  1;  . 

ablfvcU]  =  0; 
ablfevU]  =  0; 
ablratlj]  =  0; 
end; 

I* 

doj  =  1  to  7; 

if  keepO]  =  keep[j+2]  =  1  then  do; 
if  vc|j]  >  0  and  vc[j+2]  >  0  and  vc[j]  >  vcO+2] 
then  do; 

if  ((vc|j]-vc(j+2])/vc[j])  >  .15  then 
abUVc[j+2]  =  1; 
end; 

if  fev[j]  >  0  and  fevIj+2]  >  0  and  fev[j]  > 
fev[j+2]  then  do; 

if  ((fev[j]-fev|J+2])/fevU])  >  .  15  then 
ablfev|j+2]  =  1; 
end; 

if  cratiolj]  >  0  and  cratio|j+2]  >  0  and  cratiolj] 
>  cratio|j+2] 
then  do; 

if  ((cratioU]-cratio[j+2])/cratio|j])  >  .05  then 
ablrat|J+2]  =  1; 
end; 
end; 
end; 

doj  =  1  to  6; 

if  keep[j]  =  keep|j+3]  =  1  then  do; 
if  vc[j]  >  0  and  vc[j+3]  >  0  and  vc[j]  >  vc(j+3] 
then  do; 

if  ((vc[j]-vc[j+3])/vc|j])  >  .15  then 
ablfvc(j+3]  =  1; 
end; 

if  fev[j]  >  0  and  fevfj+3]  >  0  and  fevlj]  > 
fev[j+3]  then  do; 

if  ((fevU]-fevU+3])/fev[j])  >  .15  then 
ablfevO+3]  =  1; 
end; 


if  cratioU]  >  0  and  cratio(j+3]  >  0  and  cratiolj] 

>  cratio|j+3] 
then  do; 

if  ((cratio01-cratio[j+3])/cratio[j])  >  .05  then 
ablrat|j+3]  =  1; 
end; 
end; 
end; 

doj  =  1  to  5; 

if  keeplj]  =  keep|j+4]  =  1  then  do; 
if  vc[j]  >  0  and  vclj+4]  >  0  and  vc|j]  >  vc[j+4] 
then  do; 

if  ((vc0]-vc[j+4])/vc|j])  >  .15  then 
ablfvc|j+4]  =  1; 
end; 

if  fev|j]  >  0  and  fev[j+4]  >  0  and  fevlj]  > 
fev[j+4]  then  do; 

if  ((fev[j]-fev|j+4])/fev(j])  > .  15  then 
ablfev[j+4]  =  1; 
end; 

if  cratiojj]  >  0  and  cratio|j+4]  >  0  and  cratiolj] 

>  cratiolj+4] 
then  do; 

if  ((cratio]j]-cratio[j+4])/cratio|j])  >  .05  then 
ablrat]j+4]  =  1; 
end; 
end; 
end; 

doj  =  1  to  4; 

if  keeplj]  =  keep|j+5]  =  1  then  do; 
if  vc^]  >  0  and  vc|j+5]  >  0  and  vclj]  >  vc[j+5] 
then  do; 

if  ((vclj]-vc|j+5])/vc|j])  >  .15  then 
abUvc|j+5]  =  1; 
end; 

if  fevlj]  >  0  and  fevtj+5]  >  0  and  fev[j]  > 
fev|j+5]  then  do; 

if  ((fev|j]-fev[j+5])/fevlj])  >  .15  then 
ablfevtj+5]  =  1; 
end; 

if  cratiolj]  >  0  and  cratio|j+5]  >  0  and  cratiolj] 

>  cratiolj+5] 
then  do; 

if  ((cratio[j]-cratio|j+5])/cratio|j])  >  .05  then 
ablratlj+5]  =  1; 
end; 
end; 
end; 


doj  =  1  to  3; 


if  keeplj]  =  keep[j+6]  =  1  then  do; 
if  vcy]  >  0  and  vc|j+6]  >  0  and  vc|j]  >  vc|j+6] 
then  do; 

if  ((vc[j]-vc[j+6])/vc[j])  >  .15  then 
ablfvc[j-i-6]  =  1; 
end; 

if  fev|J]  >  0  and  fevlj+6]  >  0  and  fevU]  > 
fevU+6]  then  do; 

if  ((fev[j]-fevU+6])/fev[j])  >  .15  then 
ablfevD+6]  =  1; 
end; 

if  cratioU]  >  0  and  cratio|j+6]  >  0  and  cradojj] 

>  cratio[j+6] 
then  do; 

if  ((cratio[j]-cratio[j+6])/cratio|j])  >  .05  then 
ablratU+6]  =  1; 
end; 
end; 
end; 

doj  =  1  to  2; 

if  keepU]  =  keep|j+7]  =  1  then  do; 
if  vc[j]  >  0  and  vc[j+7]  >  0  and  vclJ]  >  vc[j+7] 
then  do; 

if  ((vc|j]-vc[j+7])/vc[j])  >  .15  then 
abUvc[j+7]  =  1; 
end; 

if  fevO]  >  0  and  fev|j+7]  >  0  and  fev|j]  > 
fevO+7]  then  do; 

if  ((fev[j]-fev|j+7])/fev|j])  >  .15  then 
ablfev|j+7]  =  1; 
end; 

if  cratio[j]  >  0  and  cratio[j+7]  >  0  and  cratio[j] 

>  cratio|j+7] 
then  do; 

if  ((cratio|j]-cratio|j+7])/cratio|j])  >  .05  then 
ablrat[j+7]  =  1; 
end; 
end; 
end; 

*/ 

doj  =  1  to  9; 

if  abfvc[j]  = .  then  abfvc[j]  =  0; 
if  abfev(j]  =  .  then  abfev[j]  =  0; 
if  abratQ]  =  .  then  abrat[j]  =  0; 
if  abUVc[j]  = .  then  abUvc[j]  =  0; 
if  ablfevlJ]  =  .  then  ablfev|j]  =  0; 
if  ablrat[j]  =  .  then  ablratfj]  =  0; 
end; 
run; 

data  yl987; 


seta; 

array  abfvc[*]  abfVcl-abfvc9; 
array  abfev[*]  abfevl-abfev9; 
array  abrat[*]  abratl-abrat9; 
array  ablfvc[*]  abUvcl-ablfvc9; 
array  ablfev[*]  ablfevl-ablfev9; 
array  ablrat[*]  ablratl-ablrat9; 
array  vc[*]  vc87-vc95; 
array  czone[*]  czone87-czone95; 
array  keep[*]  keepl-keep9; 
array  ab[*]  abl-ab9; 
array  sab[*]  sabl-sab9; 
array  smoke[*]  smokel-smoke9; 

j  =  i; 

if  keepD]  =  1; 
zone  =  czone87; 
countl  =  1; 
abUl  = 

abfvc[j]+abfev|j]+abrat|j]+ablfvc[j]+ablfev|j]-i-ab 

IratUl; 

if  ab|j]  >  0  then  abU]  =  1; 
else  ablj]  =  0; 

if  ab|j]  =  1  and  smokejj]  =  1  then  sab[j]  =  1; 
else  sabO]  =  0; 
run; 

proc  sununary; 
class  zone; 

var  smokel  countl  abfvcl  abfevl  abratl  abUVcl 
ablfevl  ablratl  abl  sabl; 
output  out  =  y87  siun  =  ; 
run; 

data  yl988; 
seta; 

array  abfvc[*]  abfvcl -abfvc9; 
array  abfev[*]  abfevl-abfev9; 
array  abrat[*]  abratl -abrat9; 
array  ablfVc[*]  ablfvcl-ablfvc9; 
array  ablfev[*]  ablfevl -ablfev9; 
array  ablratl*]  ablratl -ablrat9; 
array  vc[*]  vc87-vc95; 
array  czone[*]  czone87-czone95; 
array  ab[*]  abl-ab9; 
array  sab[*]  sabl-sab9; 
array  smoke[*]  smokel-smoke9; 
array  keep(*]  keepl-keep9; 
j  =  2; 

if  keepU]  =  1; 
zone  =  czone88; 
count2  =  1; 
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ab[j]  = 

abfvc[j]+abfev[j]+abrat[j]+ablfvc[j]+ablfev|j]+ab 

lrat[j]; 

if  ablj]  >  0  then  ab|j]  =  1; 
else  ab|j]  =  0; 

if  ab[j]  =  1  and  smoke[j]  =  1  then  sabU]  =  1; 
else  sab|j]  =  0; 
run; 

proc  summary; 
class  zone; 

var  smoke!  count!  abfv^c!  abfev!  abrat!  abtfvc! 
ablfev!  ablrat!  ab!  sab!; 
output  out  =  y88  sum  =  ; 
run; 

datayl989; 

seta; 

array  abfvc[*]  abfv^cl-abfvc9; 
array  abfev[*]  abfevl-abfev9; 
array  abrat[*]  abrat l-abrat9; 
array  ablfvc[*]  ablfvcl-ablfvc9; 
array  ablfev[*]  ablfevl-ablfev9; 
array  ablrat[*]  ablratl-ablrat9; 
array  vc[*]  vc87-vc95; 
array  czone[*]  czone87-czone95; 
array  ab[*]  abl-ab9; 
array  sab[*]  sabl-sab9; 
array  smoke[*]  smoke l-smoke9; 
array  keep[*]  keepl-keep9; 
j  =  3; 

if  keep[j]  =  1; 
zone  =  czone89; 
counts  =  1; 
abU]  = 

abfvcO]+abfev[j]+abrat[j]+ablfvc[j]+ablfev[j]+ab 

lrat[j]; 

if  ablJ]  >  0  then  ab[j]  =  1; 
else  ab[j]  =  0; 

if  ab[j]  =  1  and  smoke[j]  =  1  then  sablj]  =  1; 
else  sab|j]  =  0; 
run; 

proc  summary; 
class  zone; 

var  smokeS  counts  abfvcS  abfevS  abratS  ablfvcS 
ablfevS  ablratS  abS  sabS; 
output  out  =  y89  sum  =  ; 
run; 

datayl990; 

seta; 

array  abfv^c[*]  abfvcl-abfvc9; 
array  abfev[*]  abfev l-abfev9; 


array  abrat[*]  abratl-abrat9; 
array  ablfvcf*]  ablfvcl-ablfvc9; 
array  ablfev[*]  ablfevl-ablfev9; 
array  ablrat[*]  ablratl-ablrat9; 
array  vc[*]  vc87-vc95; 
array  czone[*]  czone87-czone95; 
array  ab[*]  abl-ab9; 
array  sab[*]  sabl-sab9; 
array  smoke[*]  smokel"Smoke9; 
array  keep[*]  keepl-keep9; 
j  =  4; 

if  keepU]  =  1; 
zone  =  czone90; 
count4  =  1; 
ab|j]  = 

abfvc[j]+abfev[j]+abrat|j]+ablfvcU]+ablfev|j]+ab 

lrat[j]; 

if  ab|j]  >  0  then  ab[j]  =  1; 
else  abjj]  =  0; 

if  ab|j]  =  1  and  smoke[j]  =  1  then  sabU]  =  1; 
else  sab[j]  =  0; 
run; 

proc  summary; 
class  zone; 

var  smoke4  count4  abfvc4  abfev4  abrat4  ablfvc4 
ablfev4  ablrat4  ab4  sab4; 
output  out  =  y90  sum  =  ; 
run; 

datayl991; 

seta; 

array  abfvc[*]  ablvcl-abfvc9; 
array  abfev[*]  abfev  l-abfev9; 
array  abrat(*]  abratl-abrat9; 
array  ablfvc[*]  abUvcl-ablfVc9; 
array  ablfev[*]  ablfevl-ablfev9; 
array  ablrat[*]  ablratl-ablrat9; 
array  vc[*]  vc87-vc95; 
array  czone(*]  czone87-czone95; 
array  ab[*]  abl-ab9; 
array  sab[*]  sabl-sab9; 
array  smoke[*]  smokel-smoke9; 
array  keep[*]  keepl-keep9; 
j  =  5; 

if  keep[j]  =  1; 
zone  =  czone91; 
count5  =  1; 
ablj]  = 

abfVc[j]+abfev[j]+abrat[j]+ablfvc|j]+ablfev[j]+ab 

IratUl; 

if  ab[j]  >  0  then  abU]  =  1; 
else  ab|j]  =  0; 
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if  ab[j]  =  1  and  smoke[j]  =  1  then  sab[j]  =  1; 
else  sab[j]  =  0; 
run; 

proc  summary; 
class  zone; 

var  smoke5  counts  abfv^cS  abfevS  abratS  ablfvcS 
ablfevS  ablratS  abS  sabS; 
output  out  -  y91  sum  =  ; 
run; 

datayl992; 
set  a; 

array  abf^c[*]  abfvcl-abfvc9; 
array  abfev[*]  abfevl-abfev9; 
array  abrat[*]  abratl-abrat9; 
array  ablfs^c[*]  ablfvcl-ablfvc9; 
array  ablfev[*]  ablfevl-ablfev9; 
array  ablrat[*]  ablratl-ablrat9; 
array  vc[*]  vc87-vc95; 
array  czone[*]  czone87-czone95; 
array  ab[*]  abl“ab9; 
array  sab[*]  sabl-sab9; 
array  smoke[*]  smokel-smoke9; 
array  keep[*]  keepl-keep9; 
j  =  6; 

if  keep[j]  =  1; 
zone  =  czone92; 
count6  =  1; 
abU]  = 

abfvc|j]+abfevU]+abrat[j]+ablfvc|j]+ablfevlj]+ab 

lrat[j]; 

if  abjj]  >  0  then  ab|j]  =  1; 
else  ab[j]  =  0; 

if  ab[j]  =  1  and  smoke[j]  =  1  then  sab[j]  ==  1; 
else  sab[j]  =  0; 
run; 

proc  summary; 
class  zone; 

var  smoke6  count6  abfv^c6  abfev6  abrat6  abLfvc6 
ablfev6  ablrat6  ab6  sab6; 
output  out  =  y92  sum  =  ; 
run; 

data  y 1993; 
seta; 

array  abfvc[*]  abfvcl-abfv^c9; 
array  abfev[*]  abfevl-abfev9; 
array  abrat[*]  abratl-abrat9; 
array  abltvc[*]  abtfvcl-ablfvc9; 
array  ablfev[*]  ablfevl-ablfev9; 
array  ablrat[*]  ablratl-ablrat9; 
array  vc[*]  vc87-vc95; 


array  ab[*]  abl-ab9; 
array  sab[*]  sabl-sab9; 
array  smoke[*]  smoke l-smoke9; 
array  czone[*]  czone87-czone95; 
array  keep[*]  keepl-keep9; 
j  =  7; 

ifkeeplj]  =  1; 
zone  =  czone93; 
count7  =  1; 
abUl  = 

abfvc[j]+abfev|j]+abrat[j]+ablfvc[j]+ablfev[j]+ab 

IratO]; 

if  ab[j]  >  0  then  ab[j]  =  1; 
else  abO]  =  0; 

if  ab[j]  =  1  and  smoke|j]  =  1  then  sab|j]  =  1; 
else  sabO]  =  0; 
run; 

proc  summary; 
class  zone; 

var  smoke?  count?  abtVc?  abfev?  abrat?  ablfVc? 
ablfev?  ablrat?  ab?  sab?; 
output  out  =  y93  sum  =  ; 
run; 

datayl994; 
set  a; 

array  abfvc[’'‘]  abfvcl-abfvc9; 
array  abfev[*]  abfevl-abfev9; 
array  abrat[*]  abratl-abrat9; 
array  ablfvc[*]  ablfVcl-ablfVc9; 
array  ablfev[*]  ablfevl-ablfev9; 
array  ablrat[*]  ablratl-ablrat9; 
array  vc[*]  vc8?-vc95; 
array  czone[*]  czone8?-czone95; 
array  ab[*]  abl-ab9; 
array  sab[*]  sabl-sab9; 
array  smoke[*]  smoke  l-smoke9; 
array  keep[*]  keepl-keep9; 
j  =  8; 

ifkeeplj]  =  1; 
zone  =  czone94; 
count8  =  1; 
abU]  = 

abfvc[jj+abfev|j]+abrat(j]+ablfvc|j]+ablfev[j]+ab 

IratUl; 

if  ab[j]  >  0  then  ab[j]  =  1; 
else  ab[j]  =  0; 

if  ab|j]  =  1  and  smoke[j]  =  1  then  sabO]  =  1; 
else  sablJ]  =  0; 
run; 

proc  summary; 
class  zone; 
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var  smokeS  counts  abfvcS  abfevS  abratS  ablfvcS 
ablfevS  ablratS  abS  sabS; 
output  out  =  y94  sum  =  ; 
run; 

data  yl995; 
seta; 

array  abfvc[*]  ab^cl-abfvc9; 
array  abfev[*]  abfevl-abfev9; 
array  abrat[*]  abratl-abrat9; 
array  abllvc[*]  ablfvcl-ablfvc9; 
array  ablfev[*]  ablfevl-ablfev9; 
array  ablrat[*]  ablratl-ablrat9; 
array  vc[*]  vc87-vc95; 
array  czone[*]  czone87-czone95; 
array  ab(*]  abl-ab9; 
array  sab[*]  sabl-sab9; 
array  smoke[*]  smoke l-smoke9; 
array  keep[*]  keepl-keep9; 
j  =  9; 

if  keeplj]  =  1; 
zone  =  czone95; 
count9  =  1; 
abD]  = 

abfvcU]+abfevO]+abrat|j]+abtfvc[j]+ablfev|j]+ab 

•rat[j]; 

if  ab[j]  >  0  then  abO]  =  1; 
else  ab[j]  =  0; 

if  ab[j]  =  1  and  smoke[j]  =  1  then  sab|j]  =  1; 
else  sab[j]  =  0; 
run; 

proc  summary; 
class  zone; 

var  smoke9  count9  abfvc9  abfev9  abrat9  abLfvc9 
ablfev9  ablrat9  ab9  sab9; 
output  out  =  y95  sum  =  ; 
run; 

data  dl  d2  d3  d4  d5  d6  d7  d8  d9  f; 

merge  y87  y88  y89  y90  y91  y92  y93  y94  y95; 

by  zone; 

array  count[*]  countl-count9; 
array  abfv^c[*]  abfVcl-abfv'c9; 
array  abfev[*]  abfevl~abfev9; 
array  abrat[*]  abratl-abrat9; 
array  ablfVc[*]  ablfvcl“ablfv'c9; 
array  ablfev[*]  ablfevl“ablfev9; 
array  ablrat[*]  ablratl-ablrat9; 
array  d[*]  dl-d9; 
array  ab[*]  abl-ab9; 
array  sab[*]  sabl-sab9; 
array  smoke[*]  smoke l-smoke9; 


doj  =  1  to  9; 

if  j  ==  1  and  countl  >  4  then  output  dl; 
else  if  j  ==  2  and  coimt2  >  4  then  output  d2; 
else  if  j  =  3  and  count3  >  4  then  output  d3; 
else  if  j  =  4  and  count4  >  4  then  output  d4; 
else  if  j  =  5  and  coimtS  >  4  then  output  d5; 
else  if  j  =  6  and  count6  >  4  then  output  d6; 
else  if  j  =  7  and  count7  >  4  then  output  d7; 
else  if  j  =  8  and  count8  >  4  then  output  d8; 
else  if  j  =  9  and  count9  >  4  then  output  d9; 
end; 
run; 
data  11; 
set  dl; 

format  pet  6.2  petsmk  6.2; 
pet  =  (abl  /  countl)  *  100; 
petsmk  =  (smokel  /  countl)*100; 
run; 

proc  sort;  by  zone;  run; 
proc  print; 

var  zone  countl  abl  pet  smokel  petsmk  sabl; 

title  ’1987  zones'; 

run; 

data  12; 
set  d2; 

format  pet  6.2  petsmk  6.2; 
pet  =  (ab2  /  count2)  *  100; 
petsmk  =  (smoke2  /  count2)*100; 
run; 

proc  sort;  by  zone;  run; 
proc  print; 

var  zone  count2  ab2  pet  smoke2  petsmk  sab2; 
title  '1988  zones'; 
run; 

data  13; 
set  d3; 

format  pet  6.2  petsmk  6.2; 
pet  =  (ab3  /  count3)  *  100; 
petsmk  =  (smoke3  /  count3)*100; 
run; 

proc  sort;  by  zone;  run; 
proc  print; 

var  zone  count3  ab3  pet  smoke3  petsmk  sab3; 

title  '1989  zones'; 

run; 

data  14; 
set  d4; 

format  pet  6.2  petsmk  6.2; 
pet  =  (ab4  /  count4)  *  100; 
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pctsmk  =  (smoke4  /  count4)*100; 
run; 

proc  sort;  by  zone;  run; 
proc  print; 

var  zone  count4  ab4  pet  smoke4  pctsmk  sab4; 

title  '1990  zones'; 

run; 

data  15; 
set  d5; 

format  pet  6.2  pctsmk  6.2; 
pet  =  (ab5  /  counts)  *  100; 
pctsmk  =  (smoke5  /  count5)*100; 
run; 

proc  sort;  by  zone;  run; 
proc  print; 

var  zone  counts  abS  pet  smokeS  pctsmk  sabS; 
title  '1991  zones’; 
run; 

data  16; 
set  d6; 

format  pet  6.2  pctsmk  6.2; 
pet  =  (ab6  /  count6)  *  100; 
pctsmk  =  (smoke6  /  count6)*100; 
run; 

proc  sort;  by  zone;  run; 
proc  print; 

var  zone  count6  ab6  pet  smoke6  pctsmk  sab6; 

title  '1992  zones'; 

run; 

data  17; 
set  d7; 

format  pet  6.2  pctsmk  6.2; 
pet  =  (ab7  /  count7)  *  100; 
pctsmk  =  (smoke7  /  count7)  *  100; 
run; 

proc  sort;  by  zone;  run; 
proc  print; 

var  zone  count7  ab7  pet  smoke7  pctsmk  sab7; 
title  '1993  zones'; 
run; 

data  18; 
set  d8; 

format  pet  6.2  pctsmk  6.2; 
pet  =  (ab8  /  count8)  *  100; 
pctsmk  =  (smoke8  /  count8)*100; 
run; 

proc  sort;  by  zone;  run; 
proc  print; 


var  zone  count8  ab8  pet  smoke8  pctsmk  sab8; 
title  '1994  zones'; 
run; 

data  19; 
set  d9; 

format  pet  6.2  pctsmk  6.2; 
pet  =  (ab9  /  count9)  *  100; 
pctsmk  =  (smoke9  /  count9)  *  100; 
run; 

proc  sort;  by  zone;  run; 
proc  print; 

var  zone  count9  ab9  pet  smoke9  pctsmk  sab9; 

title  '1995  zones'; 

run; 
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Appendix  V:  Work  Area  Description  and  Case  Number  Identifiers 


Section  V.  1 :  Work  Area  Description 


1  ZONE 

Description 

88th  CCSG,  SCTWO  Unit :  Telephone  Installation  and  Repair 

A1405 

88th  CCSG,  SCLR  Unit;  Metereological  Center 

A262E 

AFMC  Headquarters  Building 

A278A 

88th  CEG,  CEMIUE:  Pest  Management 

A830A-Q 

Hospital,  each  letter  identifies  unique  exposure  areas 

A862A 

Unknown 

A867A 

NAIC,  LGMM:  Operations  Maintenance 

A876C 

88th  CEG,  Grounds  Area  and  Pavement  Equipment  Area 

A878A 

Military  Golf  Course 

B18G 

Wright  Labs:  Experimental  Research  Branch 

B248A 

Armstrong  Labs:  Hazard  Assessment 

B24A 

Wright  Labs;  Aero-Diagnostics  Research 

B254A 

Wright  Labs:  Signature  Technology 

B36A-F 

88th  CEG,  Heat  Distribution  Area  B 

B40A 

74th  SGB,  Bioenvironmental  Engineering 

B433A 

NMRI,  Navy  Toxicology 

B470B 

AFIT,  ENP  unit:  Nuclear  Spectrum 

B490A,D 

Wright  Labs:  Experimental  Support  Branch  and  Avionics  Facility 

B4B1,D 

Wright  Labs;  Electro-Optics  Warfare 

B5B1,C 

ASC:  Machine  A/B  and  Machine  Repair 

B6000 

Firestations 

B620F 

Wright  Labs:  Microwave  Division 

B640B 

AFIT  /  ENP  unit:  Physics 

B652D 

Wright  Labs;  Material  and  Surface  Interaction 

B654A.B 

Wright  Labs:  Mechanics  Interactions  and  Non-Structural  Materials 

B65A 

Wright  Labs:  Fat  Frac/Rel  Grp 

B682A 

645th  LCMPS:  Library  of  Congress  and  Motion  Pictures 

B70A 

SA-ALC:  Aerospace  Fuels  Lab 

B745B,C,D 

88th  CEG,  Pavement  Equipment,  Water  Sewer  &  Gas,  Heat  Plant  Area  B 

B76A1,A2 

88th  CEG,  Firestations 

B770A 

88th  CEG,  Heat  Plant  Area  B 

B79A,B.C,E 

Armstrong  Labs:  Hazard  Assessment 

B838A 

Armstrong  Labs:  Occupational  Environmental  Vet  Medicine 

C101B,I 

Unknown 

C105A 

Unknown 

C106C 

Unknown 

C13D,F.O.P,R 

445th  LG,  Machine  Welding,  Propulsion,  Survival  Equipment,  Wheel  /  Tire 

C163A 

88th  CEG,  Fire  Stations 

C170A 

88th  CEG,  Heat  Plant  Area  C 

C174A 

74th  Med  Grp,  Orthopedic  Brace  Shop 

C19A,C 

88th  CEG,  Hazardous  Material  &  Waste,  Water  Treatment 

C2b6A,C,D,E 

ASC,  Aircraft  Modification  Division 
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ZONE 

Description 

C21C1 

Unknown 

C22A,B.C,I 

88th  CEG,  Cathodic  Protection,  CE  zone  A,B,  Project  Painters 

C255A 

Unknown 

C28B 

Unknown 

C29A1 

Unknown 

C4014A 

Unknown 

C4020A.B 

44th  LGMAF,  Fuel  Systems 

C4021E 

88th  ABW,  Age  Section 

C4024A 

Unknown 

C4035A 

88th  ABW,  Survival  Equipment  Repair 

C4066A 

Unknown 

C4067A 

Unknown 

C59A 

71st  Ordnance  Attachment:  Explosive  Ordnance 

C70B 

AFOSI:  AFOSI  Technical  Systems 

C883A 

88th  SPS:  Combat  Arms  Training 

C884A 

88th  ABW:  PMEL 

C89B 

DRMO  Area  C 

C91B1 

Unknown 

CE 

CE  Worker  Areas 

K1084 

88th  ABW,  Frame  Shop  /  Wood  Hobby  Shop 

K1240 

Heat  Plant  Kittyhawk  Area 

Section  V.2  Speadsheet  Calculations  for  Chart  4-1 


Case  # 


ZONE 

Year 

Total 

Tot  Abnm 

Pet 

B. 

a 

A1400 

89 

13 

5 

0.38 

0.1 

3.420651 

A1400 

90 

44 

12 

0.27 

0.1 

3.819144 

A1405 

90 

8 

0 

0.00 

0.1 

-0.94281 

A1405 

91 

22 

2 

0.09 

0.1 

-0.14213 

A1405 

92 

17 

0 

0.00 

0.1 

-1.37437 

A1405 

93 

18 

0 

0.00 

0.1 

-1.41421 

A1405 

94 

16 

0 

0.00 

0.1 

-1 .33333 

A262E 

90 

16 

1 

0.06 

0.1 

-0.5 

A262E 

91 

8 

1 

0.13 

0.1 

0.235702 

A278A 

89 

6 

1 

0.17 

0.1 

0.544331 

A278A 

90 

17 

6 

0.35 

0.1 

3.476344 

A278A 

91 

13 

0 

0.00 

0.1 

-1.20185 

A278A 

92 

14 

2 

0.14 

0.1 

0.534522 

A278A 

93 

13 

4 

0.31 

0.1 

2.496151 

A278A 

94 

14 

1 

0.07 

0.1 

-0.35635 

A830A 

90 

9 

1 

0.11 

0.1 

0.111111 

A830A 

91 

15 

2 

0.13 

0.1 

0.430331 

A830A 

92 

10 

1 

0.10 

0.1 

0 

A830A 

93 

19 

1 

0.05 

0.1 

-0.68825 

A830A 

94 

11 

1 

0.09 

0.1 

-0.1005 

Case  #  ZONE 


21  A830A 


22 A830C 


23  A830C 


24 A830C 


25  A830C 


26  A830C 


27A830D 


28 A830F 


29 A830F 


30 A830F 


31  A830Q 


32A830Q 


33 A830Q 


34  A862A 


35  A862A 


36  A862A 


37  A867A 


38 A867A 


39  A867A 


40  A876C 


41  A876C 


42  A876C 


43 A876C 


44  A876C 


45  A878A 


46  A878A 


47A878A 


48 A878A 


49ALL 


50  ALL 


51  ALL 


52ALL 


53ALL 


54  ALL 


55ALL 


56  B18G 


57B248A 


58B24A 


59B24A 


60B254A 


61IB36A 


B36A 


63B36A 


64B36A 


65B36A 


66B36A 


£ 

Zi 

0.1 

0.19245 

0.1 

-1.05409 

0.1 

-0.94281 

0.1 

-1.20185 

0.1 

-0.95258 

0.1 

-0.62854 

0.1 

-1 .05409 

0.1 

5.9017 

0.1 

2.933553 

0.1 

3.265986 

0.1 

-0.74536 

0.1 

0.544331 

0.1 

0.745356 

0.1 

0.745356 

0.1 

-0.8165 

0.1 

-0.8165 

0.1 

1.637846 

0.1 

0.544331 

0.1 

3.265986 

0.1 

2.236068 

0.1 

5.019011 

0.1 

1 .60591 

0.1 

0.83395 

0.1 

0.111111 

'  0.1 

0.544331 

1  0.1 

0.745356 

1  0.1 

0.745356 

t  0.1 

0.377964 

i  0.1 

7.834617 

1  0.1 

13.4624 

i  0.1 

5.885719 

0.1 

3.902887 

1  0.1 

5.003401 

1  0.1 

4.666667 

1  0.1 

4.893502 

0.1 

0.544331 

0.1 

0.235702 

0.1 

1.414214 

'  0.1 

-0.94281 

0.1 

-0.8165 

0.1 

1.382047 

0.1 

5.352644 

0.1 

2 

0.1 

2.236068 

0.1 

0.64715 

0.1 

3.666667 
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Case  #  ZONE 


67  B36A 


68  B36B 


6gB36B 


70  B36C 


71  B36F 


72B40A 


73B40A 


74B40A 


75  B40A 


76  B40A 


77B40A 


78  B433A 


79  B433A 


80B433A 


81  B433A 


82  B433A 


83  B470B 


84  B470B 


85  B470B 


86  B470B 


87B490A 


88  B490A 


89  B4gOO 


90B490D 


91  B4B1 


92B4B1 


93  B4B1 


94B4B1 


95  B4B1 


96  B4B1 


97B4D 


98  B4D 


99  B4D 


100  B4D 


101  B4D 


102  B4D 


103  B5B1 


104B5B1 


105  B5C 


106  B5C 


107B5C 


108  B5E 


109B5H 


110B5J1 


111  B5J1 


112  B5J1 


£ 

Zi 

B  0.1 

2.592725 

D  0.1 

1.290994 

1.414214 

0.1 

0.745356 

0.1 

-0.19245 

0.1 

-0.07274 

0.1 

0.261488 

0.1 

-0.65855 

0.1 

0.342997 

0.1 

-0.74536 

0.1 

-0.74536 

0.1 

-0.8165 

0.1 

-0.74536 

0.1 

-1 

0.1 

-1.44444 

0.1 

-1.17607 

0.1 

4.21637 

0.1 

1.054093 

0.1 

0 

0.1 

-0.74536 

0.1 

3.420651 

0.1 

4.924685 

0.1 

1.637846 

0.1 

0.745356 

0.1 

0.745356 

0.1 

3.444444 

0.1 

5.333333 

1  0.1 

0.333333 

'  0.1 

-0.35635 

0.1 

0.111111 

;  0.1 

2 

i  0.1 

5.851559 

0.1 

0.076472 

0.1 

-0.43033 

1  0.1 

2.434322 

1  0.1 

1.732051 

.  0.1 

1.905159 

0.1 

0.377964 

0.1 

0.430331 

1  0.1 

-1.10554 

'  0.1 

-1 

'  0.1 

3.72678 

0.1 

-0.8165 

0.1 

-0.8528 

0.1 

0.261488 

0.1 

0.486534 
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Case  #  I  ZONE 


113B5J1 


114B5J1 


115  B5J1 


116B6000 


117  B6000 


118  B6000 


119B6000 


120B6000 


121  B620F 


122B640B 


123B640B 


124  B652D 


125B654A 


126  B654B 


127B654B 


128B654B 


129B65A 


130B65A 


131  B65A 


132B65A 


133  B682A 


134B682A 


135B682A 


136  B682A 


137B682A 


138B682A 


139B70A 


140B745B 


141  B745B 


142B745B 


143B745C 


144B745D 


145  B745D 


146B745D 


147B745D 


148  B76A1 


149B76A1 


150B76A1 


151  B76A1 


152B76A1 


153B76A1 


154B76A1 


155B76A2 


156B770A 


157B770A 


158B770A 


B. 

Zi 

0.1 

-0.27217 

0.1 

-0.80013 

0.1 

3.299832 

0.1 

0.261488 

0.1 

1.60591 

0.1 

2.496151 

0.1 

0.377964 

0.1 

0.544331 

0.1 

0 

0.1 

-0.74536 

0.1 

0.745356 

0.1 

-0.74536 

0.1 

-0.8165 

0.1 

2.236068 

0.1 

-0.8165 

0.1 

0.544331 

0.1 

4.098006 

0.1 

2.370629 

0.1 

-0.5 

0.1 

0.544331 

0.1 

0.377964 

0.1 

1.222222 

0.1 

1.054093 

1  0.1 

1.054093 

1  0.1 

0.904534 

1  0.1 

-1.20185 

i  0.1 

0.64715 

1  0.1 

-0.8165 

1  0.1 

3.771236 

i  0.1 

1.905159 

0.1 

0.377964 

1  0.1 

0.745356 

1  0.1 

1.571651 

1  0.1 

2.316264 

i  0.1 

-0.88192 

:  0.1 

1.222222 

0.1 

0.777778 

0.1 

0.608581 

0.1 

0.942809 

0.1 

1.425393 

0.1 

-0.35635 

0.1 

0.7698 

0.1 

2.236068 

0.1 

2.981424 

0.1 

7.137743 

0.1 

0.342997 
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Case  #  ZONE 


159B770A 


160B770A 


161  B770A 


162  B770A 


163  B79A 


164  B79A 


165  B79A 


166  B79A 


167  B79A 


168B79B 


169  B79B 


170  B79B 


171  B79B 


172B79C 


173  B79E 


174  B838A 


175  B838A 


176  C101B 


177  Cl 01 1 


178C101I 


179 


180  C105A 


181  C106C 


182  Cl  3D 


183  C13D 


184  Cl 3D 


185C13F 


186C13F 


187C130 


188C13P 


189C13R 


190C13R 


191  C13R 


192C13R 


193  C13R 


194  C163A 


195  C163A 


196C163A 


197C163A 


198  C163A 


199C163A 


200C163A 


201  C170A 


202C170A 


203C174A 


204C19A 


c 

Zi 

0.1 

2.904382 

0.1 

1.40859 

0.1 

2.146958 

0.1 

1 .382047 

0.1 

-1 .56347 

0.1 

-1.37199 

0.1 

-0.95258 

0.1 

-1.56347 

0.1 

0.333333 

0.1 

0.235702 

0.1 

0.235702 

0.1 

-1.05409 

0.1 

-0.74536 

0.1 

-0.8165 

0.1 

-0.56592 

0.1 

-0.8165 

0.1 

-0.74536 

0.1 

-0.74536 

0.1 

-0.74536 

0.1 

-0.74536 

0.1 

-0.74536 

0.1 

1.637846 

0.1 

1.054093 

I  0.1 

1.637846 

!  0.1 

1.222222 

'  0.1 

4.157609 

1  0.1 

-0.74536 

1  0.1 

i  3.265986 

1  0.1 

-0.74536 

1  0.1 

2.236068 

1  0.1 

-0.88192 

1  0.1 

2.333333 

0.1 

0.942809 

1  0.1 

0.745356 

i  0.1 

3.162278 

i  0.1 

0.64715 

:  0.1 

4.434052 

0.1 

1.414214 

0.1 

-0.75057 

0.1 

0.273096 

0.1 

0.351364 

0.1 

3.100868 

0.1 

0.333333 

0.1 

-0.94281 

0.1 

-0.74536 

0.1 

1.728483 

60 


Case  # 


ZONE 

Year 

)C19A 

91 

C19C 

94 

C206A 

93 

C206A 

94 

C206C 

93 

C206D 

90 

C206D 

93 

C206E 

90 

C206E 

91 

C206E 

92 

C206E 

93 

C21C1 

91 

C22A 

90 

IC22A 

91 

IC22A 

92 

C22A 

93 

C22A 

94 

C22A 

95 

C22B 

90 

C22B 

91 

C22C 

90 

C22I 

90 

C22I 

91 

C22I 

92 

C22I 

93 

C22I 

94 

C255A 

90 

C28B 

90 

C28B 

91 

C28B 

92 

C28B 

93 

C28B 

94 

C29A1 

89 

C29A1 

90 

C29A1 

91 

C29A1 

92 

C29A1 

93 

C29A1 

94 

C4014A 

90 

C4020A 

89 

C4020A 

90 

C4020A 

91 

C4020A 

92 

C4020A 

93 

C4020B 

94 

C4020B 

95 

Tot  Abnm 


3 


£ 

zi 

0.1 

1.909572 

0.1 

-0.8165 

0.1 

-1.05409 

ilHKl 

-0.94281 

0.1 

0.111111 

0.1 

-0.88192 

0.1 

0.111111 

0.1 

-1 

0.1 

-0.19245 

0.1 

-0.88192 

0.1 

-0.74536 

0.1 

0.534522 

0.1 

3.01232 

0.1 

1.425393 

0.1 

0.235702 

0.1 

1.425393 

immsi 

1.571651 

0.1 

0.544331 

0.1 

2 

0.1 

0.430331 

0.1 

-0.1005 

0.1 

2.108185 

0.1 

1.732051 

0.1 

0.904534 

0.1 

4.098006 

0.1 

1.637846 

0.1 

1.290994 

0.1 

0.745356 

0.1 

-0.88192 

0.1 

0.377964 

0.1 

0.377964 

0.1 

0.377964 

0.1 

4.157609 

0.1 

-0.27735 

0.1 

0.64715 

0.1 

-0.19245 

0.1 

-0.1005 

0.1 

3.162278 

0.1 

1.905159 

0.1 

2.236068 

0.1 

4.157609 

0.1 

-0.88192 

0.1 

0.111111 

0.1 

0.235702 

0.1 

-0.74536 

0.1 

-0.74536 
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Case  # 


ZONE 

Year 

Total 

Tot  Abnm 

Pet 

£ 

Zi 

C4021E 

93 

6 

0 

0.00 

0.1 

-0.8165 

2C4021E 

94 

5 

1 

0.20 

0.1 

0.745356 

C4024A 

90 

5 

0 

0.00 

0.1 

-0.74536 

C4035A 

91 

9 

2 

0.22 

0.1 

1.222222 

C4035A 

92 

5 

2 

0.40 

0.1 

2.236068 

C4035A 

93 

6 

2 

0.33 

0.1 

1.905159 

C4066A 

91 

9 

1 

0.11 

0.1 

0.111111 

C4066A 

92 

8 

0 

0.00 

0.1 

-0.94281 

C4066A 

93 

14 

1 

0.07 

0.1 

-0.35635 

C4066A 

94 

12 

4 

0.33 

0.1 

2.694301 

C4067A 

92 

6 

0 

0.00 

0.1 

-0.8165 

C4067A 

93 

6 

2 

0.33 

0.1 

1.905159 

JC59A 

93 

6 

0 

0.00 

0.1 

-0.8165 

IC59A 

94 

11 

0 

0.00 

0.1 

-1.10554 

C59A 

95 

11 

0 

0.00 

0.1 

-1.10554 

)C70B 

90 

10 

2 

0.20 

0.1 

1.054093 

^C70B 

91 

8 

2 

0.25 

0.1 

1.414214 

SC70B 

92 

7 

1 

0.14 

0.1 

0.377964 

)  C883A 

94 

5 

3 

0.60 

0.1 

3.72678 

)  C884A 

91 

6 

1 

0.17 

0.1 

0.544331 

I  C89B 

90 

5 

0 

0.00 

0.1 

-0.74536 

!C89B 

91 

7 

0 

0.00 

0.1 

-0.88192 

IC89B 

92 

11 

0 

0.00 

0.1 

-1.10554 

iC89B 

94 

24 

1 

0.04 

0.1 

-0.95258 

;C91B1 

92 

5 

2 

0.40 

0.1 

2.236068 

i  C91B1 

93 

5 

0 

0.00 

0.1 

-0.74536 

’C91B1 

94 

5 

3 

0.60 

0.1 

3.72678 

iCE 

95 

8 

1 

0.13 

0.1 

0.235702 

IK1084 

93 

6 

2 

0.33 

0.1 

i. 9051 59 

I  K1240 

89 

29 

12 

0.41 

0.1 

5.632759 

K1240 

90 

33 

17 

0.52 

0.1 

7.949546 

:  K1240 

91 

28 

6 

0.21 

0.1 

2.015811 

K1240 

92 

30 

13 

0.43 

0.1 

6.085806 

K1240 

93 

28 

9 

0.32 

0.1 

3.905633 

K1240 

94 

28 

9 

0.32 

0.1 

3.905633 

K1240 

95 

22 

5 

0.23 

0.1 

1.989873 

Appendix  VI:  Critical  Zone  Calculations 


ZoneK1240:  Total 


‘opulation 


Year 

Total 

Abnormal 

B. 

Pet 

UCL 

1-Siama 

Case  Year 

89 

29 

12 

10% 

41% 

27% 

16% 

1 

90 

33 

17 

10% 

52% 

26% 

15% 

2 

91 

28 

6 

10% 

21% 

27% 

16% 

3 

92 

30 

13 

10% 

43% 

26% 

15% 

4 

93 

28 

9 

10% 

32% 

27% 

16% 

5 

94 

28 

9 

10% 

32% 

27% 

16% 

6 

95 

22 

5 

10% 

23% 

29% 

16% 

7 

Zone  K1240:  Nonsmoker  Only  Population 


Year 

Total 

Abnormal 

B. 

Pet 

UCL 

1-Siama 

Case  Year 

89 

8 

4 

0.1 

50% 

42% 

21% 

1 

90 

9 

4 

0.1 

44% 

40% 

20% 

2 

91 

9 

2 

0.1 

22% 

40% 

20% 

3 

92 

9 

5 

0.1 

56% 

40% 

20% 

4 

93 

8 

3 

0.1 

38% 

42% 

21% 

5 

94 

9 

4 

0.1 

44% 

40% 

20% 

6 

95 

8 

3 

0.1 

38% 

42% 

21% 

7 

Zone  B36A:  Total  Population 


Year 

Total 

Tot  Abnm 

B 

Pot 

UCL 

1-Siama 

Case  Year 

89 

21 

4 

10% 

19% 

30% 

17% 

1 

90 

35 

13 

10% 

37% 

25% 

15% 

2 

91 

16 

4 

10% 

25% 

33% 

18% 

3 

92 

20 

5 

10% 

25% 

30% 

17% 

4 

93 

13 

2 

10% 

15% 

35% 

18% 

5 

94 

16 

6 

10% 

38% 

33% 

18% 

6 

95 

8 

3 

10% 

38% 

42% 

21% 

7 

Zone  B36A:  Nonsmoker  Only  Population 


Year 

Total 

Tot  Abnm 

B 

Pot 

UCL 

1-Siama 

Case  Year 

89 

8 

0 

10% 

0% 

42% 

21% 

1 

90 

9 

4 

10% 

44% 

40% 

20% 

2 

91 

3 

1 

10% 

33% 

62% 

27% 

3 

92 

2 

1 

10% 

50% 

74% 

31% 

4 

93 

2 

0 

10% 

0% 

74% 

31% 

5 

94 

0 

0 

10% 

0% 

0% 

0% 

6 

95 

3 

0 

10% 

0% 

62% 

27% 

7 

63 


ZoneB77QA:  Total  Population 


Year 

Total 

Abnormal 

P 

Pet 

UCL 

1 -Sigma 

89 

20 

6 

10% 

30% 

30% 

17% 

1 

90 

38 

17 

10% 

45% 

25% 

15% 

2 

91 

34 

4 

10% 

12% 

25% 

15% 

3 

92 

37 

9 

10% 

24% 

25% 

15% 

4 

93 

35 

6 

10% 

17% 

25% 

15% 

5 

94 

33 

7 

10% 

21% 

26% 

15% 

6 

95 

21 

4 

10% 

19% 

30% 

17% 

7 

Zone  B770A:  Nonsmoker  Only  Population 


Year 

Total 

Abnormal 

P 

Pet 

UCL 

1 -Sigma 

Case  Year 

89 

4 

0 

10% 

0% 

55% 

25% 

1 

90 

12 

6 

10% 

50% 

36% 

19% 

2 

91 

9 

0 

10% 

0% 

40% 

20% 

3 

92 

9 

1 

10% 

11% 

40% 

20% 

4 

93 

8 

0 

10% 

0% 

42% 

21% 

5 

94 

7 

1 

10% 

14% 

44% 

21% 

6 

95 

5 

0 

10% 

0% 

50% 

23% 

7 

Zone  C22. 


:  Total  Population 


Year 

Total 

Abnormal 

P 

Pet 

UCL 

1 -Sigma 

Case  Year 

90 

10 

3 

10% 

30% 

38% 

19% 

1 

91 

12 

3 

10% 

25% 

36% 

19% 

2 

92 

11 

2 

10% 

18% 

37% 

19% 

3 

93 

14 

6 

10% 

43% 

34% 

18% 

4 

94 

7 

2 

10% 

29% 

44% 

21% 

5 

Zone  C22 


;  Nonsmoker  Only  Population 


Year 

Total 

Abnormal 

P 

Pet 

UCL 

1 -Sigma 

Case  Year 

90 

5 

2 

10% 

40% 

50% 

23% 

1 

91 

7 

2 

10% 

29% 

44% 

21% 

2 

92 

5 

0 

10% 

0% 

50% 

23% 

3 

93 

6 

2 

10% 

33% 

47% 

22% 

4 

94 

2 

0 

10% 

0% 

74% 

31% 

5 

64 


Zone  Cl 63  A:  Total  Population 


Year 

Total 

Tot  Abnm 

P 

Prt 

UCL 

1-Siama 

Case  Year 

89 

13 

2 

10% 

15% 

35% 

18% 

1 

90 

38 

12 

10% 

32% 

25% 

15% 

2 

91 

50 

8 

10% 

16% 

23% 

14% 

3 

92 

57 

4 

10% 

7% 

22% 

14% 

4 

93 

73 

8 

10% 

11% 

21% 

14% 

5 

94 

90 

10 

10% 

11% 

19% 

13% 

6 

95 

65 

14 

10% 

22% 

21% 

14% 

7 

Zone  Cl 63  A;  Nonsmoker  Only  Population 


Year 

NonSmk 

Abn  NS 

B. 

Pet 

UCL 

1-Siama 

Case  Year 

89 

8 

1 

10% 

13% 

42% 

21% 

1 

90 

16 

5 

10% 

31% 

33% 

18% 

2 

91 

24 

4 

10% 

17% 

28% 

16% 

3 

92 

24 

0 

10% 

0% 

28% 

16% 

4 

93 

33 

1 

10% 

3% 

26% 

15% 

5 

94 

43 

2 

10% 

5% 

24% 

15% 

6 

95 

28 

4 

10% 

14% 

27% 

16% 

7 

65 


Appendix  VII:  A  Proportions  Test 

The  following  test  is  based  on  the  methodology  in  Hoel’s  Introduction  to  Mathematical 
Statistics  (11:156-1 571: 


Ps  = 


514 

2454 


287 

2090 


p,  =  0.209  =  0.137  p,  -  p^  =.072 


801 

4544 


p  =  0.1763 


1 


1 


pS-pN 


^P(l  ^^(2454’^  2090^ 


0.011343 

_  Ps  ~  Pn 

^  pS-pN 

_  0.072 

~  0.011343 

Z*  =  6.34 

Zo.998  =  3.09  (1 1 :418)  and  since  Z*  >  Zo.99g  we  reject  the  Null  Hypothesis 

that  Ps  =  pN. 
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